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SYNOPSIS 


In tliic thocls, some eln.wtic inul.uoion problomo 
in two-dimenaional re^ioiia and point force .problems in 
tkcee-diDiGnsiomil rocionr. in elasticity tlxuory hitve been 
studied. Ilio tliosin is divided into two parte. Part A 
is concerned with some two-diiuexJsionai inclusion problems. 
Part B deals with x^oint force problems . 

The inclusion is a region in an elastic mcdiim 
(called tile matrix) . {CMs region (inclusion) v/ould undergo 
a. prescribed deformation in the absence of the surrounding 
material (the matrix) . But bocauae of the constraints of 
the matrix, otresses will develop everywhere in the 
material. Ihe problem is to find the elastic field in 
the inclusion and the matrix. 

Part A deals with su.oh problems. In the works of 
the previous authors, the inclusions wore supposed to be 
present in an infinite or semi- infinite medium. Ihe 
methods of solutions of such problems were based upon 
’point force’ technique or energy considerations. If 
the matrix is finite these methods do not succeed 
because point force requires an infinite or semi- Infinite' 
medium and energy considerations require considerable 
guesBing of the equilibrium interface. 



Practically no problems woi’o earlier solved for 
f ini, to rogiono. Fortunately, bowevor, Mioa,i the matrix 
is finite, tlie tboory of Etlbort problem and complex 
variable metbod can bo applied. As far as is known to 
the author, tlxi-s method has boon applied for the first 
t.lmc‘ in tills thesis after suitable adaptation. Bxact 
analytical solutions have been obtained for some problems. 
IhosG probloras are (i) an occeiitrio circular inclusion in 
a circular region, (il) an olUptic inclusion in a circular 
region and (ill) rectangular inclasion in it circular region. 
Phe solutions can be adopted to the case v/iien the Poisson 
ratios of the inclusion and matrix are different but shear 
moduli the same. It niay be noted that the results obtained 
by other authors for infinite medium can be obtained as 
particular cases of the solutions obtained in this thesis 
by suitable limiting processes. Ihis has been indicated in 
this thesis. 

As stated above, the point force techaique has proved 
to be quite useful when inclusions are present in an infinite 
medium. However, it may be mentioned that the point force 
is not only useful for inclusion problems but even otherwise 
it has found applications in material testing, soil meolmnios 
and engineering structures. Ihus the results for point foroe 
are interesting in themselves. Hence in the. Part B of the 



thesis , we have solved some point force (also called 
conoentratod force) problems , 

Ihe first problem is that of a cone under an axial 
point force acting at a point on the axis other than the 
vertex of the cone. It may be mentioned that the problem 
of axial point force acting at the vertex of a cone was 
solved by J,H. Michell in 1900 and since then no progress 
was made. Iho results of our problem involve infinite 
integrals which are convorgent , lo have the feeling for 
the stresses in the cone for this case, the integrals were 
evaluated numerically for some cases and have been given 
in the thesis at suitable places. When the point force 
acts very near the vertex, it is obvious on pliysical 
grounds that the stresses cannot differ significantly from 
the case when the point force acts at the vertex. Numerical 
results confirm this escpectation and are reported in the 
thesis . 

Ihe second problem is that of an infinite cylinder. 
Here again the axial point force acts at a point on the 
axis of the cylinder, which has been chosen as the origin. 
!Che method is essentially the same as that employed for 
the cone. Ihe results involve infinite integrals again. 
Ihese are evaluated numerically and presented in the thesis. 
It may be noted ttiat when the radius of the cylinder is 
very large , we deduce the case of a point force in an 
infinite medium. 



The next prolalem is tlK.it of a point force in 
a wodge . Ihe point force acts at any point (not 
necessarily on the axis of the v/edge) and is inclined 
to the axis of the wedge at any angle . fhe results are 
tabulated for some particular cases. 

In all the cases of point force, an auxiliary 
problem when tractions are prescribed on the surfaces 
in the case of cone and cylinder and on tho faces in the 
case of vrodge, is also solved. 

Lastly the problem of a circular inclusion in a 
wedge is solved. Results are given for a symmetrically 
situated inclusion. Ihe problem of a circular inclusion 
in a half plane and an infinite medium can be deduced as 
particular cases. 
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lirfRODuai'ioii 


Au in tiiQ (j^'-jiiopLiic , thiu tlie^io doola 

witli aoLio iucliujion probliniiu in tv/o--diiflo.uuion!'; and point 


forco problomy in tlirGO-dimoiioiono 


in tl'io clo,.-'. 


iciil 


irjotropic olanticitir' tboory . 


Tlio doj'inition:! and tlio 


pro'bioiiiL; aro aloo 


otated in t.Uo Liynupoio. 


Tlio otudy oi' iuu-l.uoion problomo wan t'irot 
initiated by Prciilcol ((25)) in coimoction v/itli liic' 

'kinutic tlioary of liquids and by Mott and Mabarro ((26)) 
in contieotion with their theory of precipitation hardonind 
of alloye . litirther pro^rooo was laade only after iSohelby’c 
papers ((11, 14)) dealing v/ith the elastic field of 
ellipsoidal inclusion* 

Since then rapid progroso has been made. Reference 
may bo made to the classical v/ork of Jaswon and Bhargava 
((15)), where for two-dimensional problems the point force 
was coupled with complex variable technique to obtain 
esqplioit solutions for elliptic inclusion problem. Ih© 
technique could be further applied to anistropic bodies ^ 
and this was done by Willis in ((£?)) and by Bhargava and 
Kapoor ((28, 9)) to obtain tho solution for the triangular 
and rectangular inclusions in an infinite medium. I'he 
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latter ;'.utiioi.'e aluu tjolvcd tjoiau iiioi’e intoi’cat-Uiij: problems 
ron.o.tiiiu to iiiuluiJioiiu a.u iialf jalnne and tlie probloitm 
vdioro tile incluuioii ixitcractud \Yitli anothor iiiclnyion or 
an inlioraobcnoity or a cavity in an othorwioo infinite 
modiura, SliQ.nua ((29)) considered tlio circvilar incluaion 
in an infinite atrip . 

In 1961, Bliaroava ((lo)) ahowed tint tlio principle 
of minimtun atra.in energy can o,lso be aiUiccaafuHy employ od 
and ia more useful v/lien in place of incleiaion tliuro is an 
inliomoccnoity . Ihia idea was used by Blia.rtjava, and 
Hadhalcrialma. to obtain aonio G::plicit aolutionc for isotropic 
((12)) ?.aid aniotropic ((13)) elliptic inlioitiooonoitios in 
an infinite medium. Bundura and Mura ((30)) iiavo investigate' 
tbo beliaviour of an edge dislocation present in tlie matrix 
witli a circular inclusion. Bundurs and Sendeckyj ((3l)) 
have studied a aimilar problem in v^hlch the dislocation 
is present inside the circular inclusion. Applications to 
some of the results t'o iron nitride precipitates in iron 
were made by Bhargava and Mclean ((32)). 

Knoxes ((33)) derived on etiuation for the strains 
of an arbitrary elastic field in an infinite matrix 
perturbed by several inclusions and solved the problem 
when the shear moduli of Inhoraogeneities and matrix are 
identical. Some other recent contributions .in tMs field 



of Qtudy liuvo iDoon coiicorned vdtli ULiinj Yfiriationnl itiotliodc 
to dorivu bomuLa for tlio aiji.iro{so/t;o iiioduli of inu'lttpho.cod 
matoriJila Ikaviri^i arbitrary phaoe geometry. ((b4)) 

octimatod tiio ovoralU. moduli of an arbitrary fibre conipocite 
v/itli troncvoroGly iootropic plia.Guo and aloo tlio macroocqpio 


elaotic moduli of tv/o phase comxjositos ((oS)). Jiudiansly 
((3^)) analysis for the dotoriaina.tion of tiio olaotic 

moduli of a compoaito material , Ihe bounds for olantic 
fiioduli. of solid composite matox’ivals wore {jivon by Walpole 
((37, 38)) by employing; e:d; remum principles. 

In o.i:i the problems Biontioned above, the matrix 
was tfdceu as infinite or somi-ixrCinite . As far as is known 
to the author, inclusion probli-'ras v/hen tho matri:-: is finite, 
has been solved only in this thesis. Such pro blows are more 
roalistic and teohnoloeically important. 

As mentioned earlier, the exqxlicit solutions vrore 
obtained either by point force teclmiipue or energy 
CO ns id orations. The point force requires an InfiirLtc or 
semi-infinite modia and energy considorations require 
considerable guessing of the equilibrium interface, both 
of which are not easily possible if tho matrix is finite. 
However, the theory of Hilbert problem coupled vYith complex 
variable can bcb'.qn)licd to solve some' inclusion problems 
which are reported in Part A of the thesis . 
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Part A consist t3 of four diaptero . In the first 
chapter some results of the complex varialDle formulation 
of two-dimensional elastostatic problems which are needed 
in subsequent chapters are given. Also theory of Hilbert 
problem and the method of application of this theory to 
inclusion problems is briefly recoiinted in this chapter. 

In chapter II, the problem of eccentric circular 
inclusion in a circular ring is solved. Exact analytical 
solution has been obtained. The results of concentric 
circulo.r inclusion in a circular region, circular inclusion 
in an infinite medium and circular incl\ioion in a serai~ 
infinite medium are obtained as particular cases of our 
problem. Variation of stresses at the eq.uilibriUDi boundary 
both for inclusion and matrix is shovm for some caaos in the 
form of graphs v/hich are given in Appendix to chapter II. 

Olmipter III deals with the problem of concentric 
elliptic inclusion in a circular region. The results of 
elliptic inclusion in an infinite) medium and the results of 
concentric circular inclusion may be derived as particular 
cases of our problem, lines of maximum shearing stress 
have been drawn for matrix and their grapte are given in 
Appendix to chapter 111. 

The problem of rectangular inclusion in a circular 
region is solved in chapter IV. An exact analytical 



i.Jo'Uitioii has hoon i^lYon. Sq^uaro iiioliisionc in a 
circular aiitl infinite ro^’icn and rectaiich-lar incltioion 
in an infinite ro£i'ion arc obteiinod as particular cases. 
For matri:c lines of iiiaximmi uhearin^j strof.ui have beGn 
drawn and are QiYon in the form of Appendix 

to chapter IV. It is intorestinq; to note that the lines 
of maximum shcariruj stress omancatin:'; from the hoimdariGS 
of slondor ci;Lipscc resemble to some extent v/ith those 
emanatinii' from the boundaries of roctanjlos with small 
width. 

As stated earlier point force tocimiquie has 
been widely used by some authors in solving; inclusion 
problems. Ihe results for point force are thus useful 
for inclusion problems. But elastic field duo to a 
point force is important by itself also and has useful 
applications in soil mechanics, material testinc and 
orji{5inc':rini; structures etc. 

The solution of the problem of a point force 
at any point in a threc-dimonsional infinite elastic 
medium was given by Lord Kelvin ((S9)) and is described 
in love’s book ((16)). ^Dhis solution has been extended 
to the case of a semi-infinite elastic body under a 
point force by Bean, Parsons, and Sneddon ((diO)) and is 
readily available in ((l7)). Michell ((41)) gave the 
solution of the problem of axial point force acting at 
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the vertex of the cone. ®he complex potent tala corresponding 
to a point force in a two-dimensional infinito elastic 
plane are given in G-reen and Zerna ((l)) and by Bhargava 
and Kapoor for a point force at the boundary of a circular 
inhomogeneity in an infinite medium ((42)) and also in a 
semi-infinite medium having a circiilar cavity ((43)). 

In Part B of the thesis some point force problems 
in three-dimensional and two-dimensional regions have been 
considered. 

The technique which is employed in the Part B 
of the thesis may be very briefly described as follows t 
Consider an infinite clastic region in which a point force 
is acting, we shall get a system of stresses everywhere which 
is being called the first stress system. Imagine the 
surface (in our cases, a cone or a cylinder or a wedge) 
within which the elastic field is to be considered. On. the 
surface of the body, tractions would be acting. Nullify 
the tractions by applying equal and opposite tractions and 
get a second system of stresses. Superpose the two stress 
systems* Ihis gives the required stresses witliin the body. 

As the problems, that we shall deal with are 
axially symmetric, the theory and some results are briefly 
described in chapter V. (Details may be found in ((16)) and 
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((l-''0) )• ^31 tliin chapter tJic following rooulto are gi-vcii: 

(i) llio 0 i|uutinnrj of oiiuilibrima in ciO-inirical polar 
coorlinates , the oolntion of which dopeiitlo upon two 
functiono, (ii) HJhe diapluciMGiit compononto duo to a point 
force in threG-dimonoional region, (iii) lEhc [jtreofjoa in a 
t\vo~diiacncional region uiido'r a point force. Ihioy are neodod 
in the oubueciuent cliai)toro. 

In cliapter VI, v/o consider tho caoc of an a;:ial 
]Doii)t force actin^^ at a point on the arJ-o of the ooni-inf Inito 
claotic cone. It may be noted that the prcvioiia solution by 
Wicholl ((41)) doalo with the proilera of a:cial point force 
acting at tlio vertex of the cone . Ac otatod earlier, firfjt 
we find the stresaes in a infinite medium duo to the point 
force. Y/g got the firot ayotem of ctroaceo. Now in this 
infinite media, imagine a cone (origin iu the vertex of the 
cono and axis of tho cone coincidoo with the direction of 
the point force acting at the point (0, 0, d) ) and find the 
tractions on its surface. We apply the tractions on the 
surface of tho cono equal and opposite to that obtained 
from first stress system. Iliis gives rise to another 
stress field in the cone which we call the second system 
of stresses. O’hese two stress systems are superposed 
and we obtain the stresses in the cone due to a point 
force with zero tractions on its surface. Ihc, results 
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for 'bho firot and aecoiul syatem of stresoeo are given 
in tliis chapter. 

In cliax^tcr VII, the mathematioG employed for 
finding the second oystem of otreGoes io given, as 
it ia not trivial. This involves the reoalts of 
chapter V . First results of cliapter V are converted 
into spherical polar coordinates . D3he two functions 
mentioned in chaxotor V satisfy oexuations (ll2) p. 
given in chapter V and when Mellin transform of those 
equations, is taken, the solution of the transformed 
e.xuations iivolves two unlcnown constants ?/hich are 
determined v/ith the help of two houndary conditions. 
Substituting the values of the corstonts in the stresses, 
we get the secorti system of stresses. Ihe 3tr0S,';eu are 
given in terms of infinite integrals. 

In chapter VIII, some graphs and tables for 
some interesting values of d (the distance of the 
point of application of the force from the vertex of 
the cone) are given. Ihis is necessary because the 
results of the second system of stresses involve 
infinite integrals and a method is needed to evaluate 
the stresses numerically. As check on the numerical 
work, one of the values of d was taken to be 0.01 
in non-dimensional form so that the point force acts 
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Tory tlio Yortex of tlie cone . It io otvious on 

pliysioal ti'rounds that the atreijyerj caimot differ 
significantly from the case when the point force acta 
at the vertex. Huuierical results confiria this 
e.-xpectation. Is the results involve infinite integrals, 
it is necessary to consider the convergence of the 
integrals which lias been done in this chapter. 

In chapter IX, the problem for an axial force 
in an infinite elastic cylinder is considered, Ihe 
method is essentially the same as described in chapter "VI. 
Iherefore only the results for the first and second 
stress systems are given, Ihe second system of stresses 
involve infinite integrals. The convergence of the 
infinite integrals is discussed. 

In Appendix to this chapter, some graphs are 
given for the stress field for various sizes of tlie 
cylinder. 

In chapter X, we deal with the point force 
problem in a wedge in two-dimensional elasticity theoxy . 
But this force acts at any point and in. any direction 
in the wedge. Thus we have considered a point force 
(Xq, y^) acting at (h, h), , The method of solution 
is the same as in chapter YI. In Appendix to this 
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chajjter, the convergence of the integrala io diociiosod 
and also some graphs and tables are given. Results 
obtained by the techniiiues of the cliaptor VIII are 
compared with the Icnown solutioais in a partictilar case 
and are in very good agreement. 

• lastly in chapter XI, we consider the problem 
of circular inclusion in an infinite elastic wedge. Ihe 
results are applicable for circular inclusion anywhere 
within the wedge. Results of circular inclusion in a 
half plane may be obtained as a particular case. 

The v/orlc in chapters II, HI* IV, VI, VII, VIII 
and IX is based upon the following papers some of which 
are published and others are under publication. 

1 , Eccentric circular inclusion in a circular region 
(Published in the Journal of Physical Society of Japan, 

Vol. 25, ho. 3, 1968). 

2 . Elliptic inclusion in a circular region (Under 
publication in Proceedings Natioiml Institute of 
Sciences, India). 

3 , Rectangular inclusion in a circular inclusion (Under 
publication in Bulletin de 1’ academie Polonaise des 
Sciences, Poland). 



4. 


0 


ill rni I'lairbic 
;it a point on t^io a'-ui 


(iono ''liio to an aainj. 
(Putli-oliiul in Ac’ba 


.i’orco 


i.ioojiojiica, Aiir.tria, Vol.O, 110,4, lU 
Strcrjtioti in an infinite- olaotic cyliU' 1 oa? - 4iio to a 
point foaxe valony tiio ardr. . (Pnliliohod in 
I.iito'rnationn.1 Journal of Enyinotjrin;'; Stdonco, 'iT,S 


Vol.7, 1 . 10 . 7 , Id 00). 
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CHAPTER I 

COMPLEX VARIABLE POm/rOLiTION MB THE HILBERT PROBLEM 

This chapter de£?cribeg the formulation of two- 
dimensional elastostatic problems by the complex variable 
approach and the problem of linear relationship (the 
Hilbert problem). This is done to medce the thesis self- 
contained . 

The method of solution of two-dimensional problems 
in infinitesimal theory of elasticity and the history of 
development have been exhaustively dealt in the books of 
Green and 2erna ((l)), Muskhelishvili ((E)) » Sokolnikoff 
((3)) and Timoshenko and Goodier ((d)) . 

In the two-dimensional elastostatic problems of 
plane strain and generalised plana stress, the equations 
of equilibrium in the absence of body forces are 
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9 P 

XX ^ 

8 P 

~ 0 

^ u j 


9x 

0y 


a p 

7^ + 

8 P 

— yy = 0 , 

(1) 

5x 

8y 


where Pj^^ ( 

i»3 « 3c,y) are the stress 

components and x,y 


is tlxe Cartesian coordinate system. 


The compatibility condition which is satisfied in 
the region R under consideration is 


V^(P + P ) « 0 . (2) 

^ ' XX yy ' 

The conditions on the boundary 0 of R are , 

^xx. cos(n,x) + Py^ coB(n,y) := 

( 3 ) 

and P^ cos(n,x) + P^y cos(n,y) =* P^ , 

where n denotes the normal to the surface element at 
(x,y). 

As noted by G. B. Airy, the eq.uations of equilibrium in 
(1) are satisfied by substituting 


p „ afs. 


P W * 
xy 


9 % 




anu ^ o * 

yy 


( 4 ) 


where B « TT(x,y) is a real function. Substituting 
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in (fi), loads to the biharMonic equation 

= 0 . ( 5 ) 

Ihn solution o.C (6) in terms of two analytic functions 
^( 25 ) aaid'Y- ( 2 ) can be written as ((3)) 


.^U - z ^ 'Y-l-?; f (s) 

where 0 =: x + iy, and bar stands for complex conjugate. 
From (4) and (6) 


■p 



afu 

_ i ^ , 


^xx 

“ ^xy 


8 

dyax 


F 

- i P , 


»aaj-.--;r- 

i 9% 

•mt ' • 

(?) 

yy 

xy 


d x^ 

dydx 


Also, 






a. 

i&TJ 

(|}(a 

) -f 0 

^)*(z) + M^(z) » 

(8) 

9x 

9y 





where dash means differentiation with respect to the 
argument throughout the thesis and vj'fe) = . 

Finally, the stresses and displacements in 
Cartesian coordinate system are given by 


^XX ^yy ” 

^yy - ^xx + V “ ® 


( 9 ) 


and5,jJ-(u tiv) « K |)(3) ~ z (|)Us) - ’ 


(10) 


Re stands for the real part of a complex quantity, 


where 
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K ™ 0 — 1 - "v) .fo.r ctr-'i-in -'viid 'iC = for (>'cn(irolLi:.3GcI 

pltmo '0 'blic Poioo.on ratio. 

I’lio analytic functions (|)(s) and vjj (&;) are 
doterminod unioucly for a {jiven otate of otreoo, if tlie 
ori(;in of cocjrdiriatos io talcen within tho rOiiiion tuider 
conL^idoration a.nd the followinc; conditions are imposed 


\>{Q) ^ 0 , I(<1)'(0)} =0 and v|j (o) = o, (ll) 

whore I stands for the imajinary p^^rt of a complex q.iiantity . 
Por the first "boundary value proilom i.e. when the tractions 
arc preGcri"bod on the "bounda-ry 0 of R, the "boundary 
condition is 

(j)(a) + z iS)’(!3) + H^Cz) = f.(s) + i fo(£3) + Qonstant on 0, 

( 12 ) 

s 

where f^(s) + i %('^) - ^ f (^nx ^ 

Oonforraal trans format ions have been utilized effootively 
in tv/o-dimcnoional elasticity problems . A brief account 
of the conformal mai)pin{i and fhe methods of determining 
the analytic functions ^(z) and (z) are given below. 

Suppose the region R (finite or infinite) whose 
boundary 0 has continuously changing curvature is shflply 
connected and there is an analytic function 

-L ^ (13) 

which maps the region R conformally on the unit circle 
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We disting-uish a point on the ’boundary of (ft si by 
writing 5 = <»- , 

Denoting 


<1^1 (- 5 )= and c|>'(7:) ;: 


( 14 ) 




and noting that s' r on the boundary of the unit cirole 
equation (l2) may be written as 


c|) c«-) -f 




4* ^ C«^) -4- MJ 0»') 




(15) 


where H(r') « + and on the boundary of the 

unit circle . 

The expreoeions of stresses are given by the formulas 


^ Ret %CJ,3 


( 16 ) 


' V v- 


p +4t p 


A 


P "VU- 'cJ ) 


lu-Ul §(i)+ io-'cs)^«)] 


d’i’) 


where ^ c-s ) = and f 

’Vo- 'ts ; "UA '«) 

For a finite simply connected region R> the functions 

(fjCJ) and nay be represented in the power series 


>!»# 


, 4 >,«) - E H’,W)= 2: in R. (l8) 

■MSo 


Both H(-) and in (16) ©an be ©spanded by the 

’'VO' 4 .'“ ) 

oomplex Courier series as 
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_ 5" and = Z , (19) 

-U>V) -(X5 

By substituting those aeries expansions in (18) and ( 19 ) 
in (l5), and on comparine the coefficient of different 
powers of a- on both sides of the equation, the 
undetermined coefficients and b^ are obtained . 

In order to obtain unique solution in the first 
boundary value problem X <{ assigned and 

cja, Co) = 0 • 

The analytic function^ 4>C-5) and ^^Li) can also be 
determined with the help of integrodifferentlal equations* 
Multiplying both sides of (IB) by -L , where 

Jttl 

and integrating over -V the boundary of the 
unit circle, it becomes 


J_ ( 




i 


where 


X. r 

4!itC Y a*> <»“— 

/\(4] = JL. f 

JZTTi. •'Y 


/\(f) (so) 

V -"--s 


i 

By Harnack's theorem, equation (iS) and (20) are equivalent* 
The analyticity of and reduces (EO) to 

the desired integrodifferentlal equation 


+ -±X . ( 81 ) 

Th© value of the unknown constant v^fo) may be determined 
by imposing the condition ())(,o) = o. 

may be determined by multiplying the conjugate 



and integrata.ng over W . Ilie 
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of (IB) by-L ^ 

5!Tft 

required equation io 


vlil-O =, JL f 4- 

' SttC Jy 



^TTf 


OUC^) 





TMs above formulation of the two-dimensional elastostatic 
problems in terms of the analytic functions ^(a) and (g) 
shall suffice for the purpose of this thesis. 

theory of Hilbert problem is now stated below in 
brief. This theory is described in details in ((l)). I'or 
defining the Hilbert problem, we begin with the definition 
of a line. 

A line is defined as the union of a finite number 
of simple, smooth, non-intersecting arcs and contours in 
the complex z plane. The line is denoted by I. 

It is assumed that each arc or contour which is a 
component of 1 has a definite positive direction. The 
ends of the arcs, if such exist, form part of I and are 
called ends of the line L, These ends are denoted by 

Draw a circle of sufficiently small radius about any 
point t which does not coincide with one of the ends of D. 
The two parts into which the circle is divided by the line 
I are called the left and right neighbourhoods of the point 
t -on D according as they lie on the left or right of the 
line on looking in the positive direction of D. In a 
similar manner the "left and right neighbourhoods of any 
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part of the line I may be defined, flies e left and right 
neighbourhoods arc distinguished by (-f) and (-) respectiTely . 

A function is continuouc from the left or right 

at a point h on L , if FCa) tends to a limit as z tends 
to t along any path remaining however on the left or rig^it 
of L. fho limiting values are denoted by 

I'‘^(t) or P"(t) .. 

Let s’ be the z plane cut along L and F(z) be some 
function, given in S‘ but not on L and satisfying the 
following conditions 

1. fhe function F(z) is holomorphic everywhere in S’, 

2. i’(a) is continuous from the loft and from t^ie right 

at all points of L, other than the ends * 

3. Hear the ends 

A 

|p( 2 )| < , 0 s< 1< 1 , (23) 

|z~c| 

where c is any one of the ends ^ ^ constants, 

fhen such a function I?(a) is called sectionally holomorphic 
in the entire z-.plane or simply sectionally holomorphic, 
fh© line L is called the line of discontinuity of fCz), 

How the Hilbert problem is to find the sectionally 
holomorphic function I’(z) with the line of discontinuity L, 
the boundary values of which from the left and from the 
right, satisfy the condition 





= G(t) F'Ct) + f(t) on L , (24) 

orcept at tlie ondo * G-(t) and f(t) are the f-unctiona 

{ji’ven on L and Gr(t) ^ 0 evorywliore on I* GrCt) and 
f(t) oatiofy tlio conditions in (23). 

The oimplest case of the prohlom oocura when 
G-(t) = 1. Tho aiuatlon (24), novf becomes 

3i'‘*'(t) - ^“(t) = f(t) on I. . (25) 

The Golirbion of (25), may bo written as a Oaiioliy intcd'^al 


= 


1 c fCtJa-t 

27 ci h 


(26) 


It may be seen from (26) that i\(^(z) vanishes at infinity 
and satisfies (25), 

For finite regions, the solution of (25) may be 
given v/ith the help of the f ollov/ing result * 

If FCz) is oectionally holomorphio in some 
finite region Sq , then this function F(.'3 ) may always 
be represented in the form of the sum of a function Fq(z) 

holomoxphic in S and a Cauchy integral F^(z) si,/ 

Thus F(z) « ^ 0 ^^^ ■*' explicitly 

ti* ( n\ ^ 1 r f , (t , ^,-d;^. , ( 2 * 7 ) 

i a + 2%il t~Z 

li is the line of discontinuity and f(t) is given by (26). 


!Dlic expresaion in (2?) hold true eYerywhere in S^, except 
at the points of L . 

Che theory given above may be used to solve some 
two-dimensional elastostatic problems as follows. Let S 

o 

be the two-dimensional region under consideration. Che 
outer boundary of Sq is denoted by which for the 
purpose of the subseq,uent chapters, may be assumed as 
simple closed contour. Che discontinuities in the 
displacement components for a passage tlirougii L are given 
by 

u'^(t) - u’“(t) = ’ 

( 28 ) 

v'^(t) - v”(t) = ^ » 

where g^(t) and gg(t) are known functions and (+) and 
(-) superscripts are the same as defined previously. 

If and are Cartesian components of the 
external stress vector acting on the boundary , and 
n is the outward normal, than the following boundary 
conditions are satisfied. 




(){t) + t F(tT + PTC = f(t) on 


( 29 ) 
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+ g^g(t) onL,(3l) 
s 

where f(t) = I f (X^ + i Y ) ds on I (32) 

and g(t) « g^(t) + i gg(t) on 1 (33) 

are known functioi^^ and ^{z) and ^ (z) are the same as in 
(lO) and (12). 

From (29), (30) and (3l) , we obtain 



.+/ N K-/ ^ 

(ir(t) _ (t) = 

on 1 

(34) 


K + 1 



and 

^ 2/th(t) 

'j^(t) — ^^(t) S3 -H,*---.— 

on X, 

(56) 


K + 1 



where 

h(t) B - g(t) - t g'(t) , 

g-ct)= Ifhi . 

(36) 

Ihe solutions of (34) and (36) as 

given in (2?) may now 

be 

written as 




t)(e;) B + ^)#(z) 


(37) 

and 

^ (z) « %iz) + %(z) 


(38) 

where 

(|)q( 0 ) and ^q(z) are functions holomorphic In 3^ 

and 


/t f, 

<P*(z) sa ^ J 

ii;i(K+l) Xi 


C»9) 


t - 0 




a 


r SiSJa . (40) 

7i:i(K+l) L t-a 

It may be noted that the fmaotions |)(zi) and '|»( 2 ;) replace 
P(z) in (2?), and similarly (l)^(z) , ^’q(z) replace 
and it)*(z) and replace P*Cz). 

The integrals in (39) and (40 ) have tv/o values 
depending upon whether z lies inside I or outside h. 
Substituting (39) and (40) in the boundary condition (29) » 
we obtain 

Mt) + -t == ^0 

O 

where ~ fCi*) ~ ■“ ^ 

and f^(t) is a known function on . 

The functions (z) and (z) may be found out with the 
help of series expatision method or integrodifferential 

equations as described earlier* 

The results of this chapter are applied to a few 
technically important problems in the following three chapters. 


and ^*(z) 
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GHiPTER II 

EGOEMIHIC CIRGTJXAR INCLUSION 11^ A GIROIILM REGION 

Ihe problem of eccen-feriG circular inclusion in 
a circular region is considered in this chapter. She 
problem may be stated as follows : 

A circular ring of outer radius R has an eccentric 
circular hole. Ihe outer boun<iary of the oircxilar ring 
is denoted by and its inner boundary by h . Ihe 

centre of is denoted by 0* and that of Ii by 0. 

She radiu.B of Ii is r’ and c is the distance between 0 
and 0* (Figure 1 p.4s ). This medium of circular ring 
is called matrix. If an elastic solid called 'inclusion’ 
of dimensions sli^tly bigger than those of the hole, but 
remaining within the limits of proportional elasticity is 
embedded in the matrix, then because of the misfit in siae, 
stresses would defelop both in the matrix and in the 
inclusion . 
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©lo problem caii be viaualiKed ao follov/o i 

Oonj.iidor a circular rc^ilon (called Incliujion), 
occoiitrically situated in a circular rinc (called matrix), 
©lie inclucion luxdcrii’oco a cpontaiieoiiG atreas free 
transformation v/liich in the absence of the matrix is 
precoribod. But because of constraints of the matrix, the 
ctroiJijoG and atrains appear both in the matrix and inclusion. 
©iG problem io to evaluate this elastic field everywhere. 

It may be seen that the mathematical solution io the some 
for the tv/o problema stated in tliis and provious pr. rafjraplip . 

©10 incluoion in the abcence of the surro-undinL<’; 
material undor^^ooD a displacorient whoae components are 
oharacterieed by (g^x + 6^y, G^y + 6gX) with respect to the 
origin at 0 i.e. the centre of B . 

Bet (u"^, v”*") be the displacement components of the 
inclusion and (u“, v“) be those of the matrix. 

At the equilibrium interface, if the displaoeraent 
components are denoted by (uj , v^ ) and (u^ » v^ ), then 

u J - u^jj = - (g^3c + 6^y) == g^(t) , 

( 42 ) 

« v^ « - (Ggy + 6gx) = gg(t) . 

It may be remarked that (<’ < 5 , the displacement at 
the boundary of the inclusion is measured from its 
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natural state; the state it would have been in the 
absence of tlio uiatri^c; while ( \C , ) ; 


tho 


displacoKiont at the boundary of mo.trix mofi.Burod from its 
natural state, v/hich is the state vdien tho inner hole is 
circular. 


fho method of Go3.ution of tliio problem is bohod 
upon Hilbert's thoo3?em for fiaiite roGions wJiioh Jias been 
stated in chqpter I . 

Usint-; complo’v variable approach tmd Hilbort theorem, 
tho problem is reduced to tho determination of tv/o 
soctionally holomorphic functions and 'jw/(a} , having 

the line o.i discontinuity L. Each ij){ 2 ;) and vj/ (s) may be 
expressed as the sura of a function analytic in the whole 
region i.e. l 2 ~cl ^ R and Cauchy integrals evaluated 
along the circle L , Pirst g(t) and h(t) are evaluated 
from (33) and (36) and then and from (39) 

and (4o) . 

'file equation of L is jz) = 7 -^ * ^or a point on 
L, a (z =s X + iy) is denoted by t. Substituting the 
values of ^(t) and gg(t) from (42) in (33) and (36), 
v/e get 

2 

g(t) = -{Cej^+ eg) - 5g)} I - (e^-eg+KSj^+fig)} U , 

frrf U 


# • 


. (43) 
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Ji(t) = (ej^+ eg) t ~ - 

r3' 

-^e^- eg+iC6i + 5g)} r-«J_. ( 44 ) 

2t'^ 

Theoe values ol' fi(t) and li(t) arc siibotituterl in (39) and 
(40) and alter oval^uatlnc contour integrals, we got 







+ Gg - i(&4- 6g)| 

(K + l) 

/l ^Gi-Gg ~i(&^-fr6g)] 2 

(K + 1) 


( 45 ) 


/t r* >i-i(e>^+ 6g)j ^ ^ ^ 

(K + l) ’ ^ 

( 46 ) 

g/lr*^(Gj^te 2 ) ^ /I r* ' ^ 

(IC +1) z (K +1) 


where the su'bscrii)t i has been used lor the inclusion 
and subscript m for the matrix. 

fhe boundary condition satisfied on the outer 
boundary is 


« f^^(t) 


(t)0(t) + t 1)^(1) + 
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vYjicro f^Ct) = fC'fe) “* ^ '4^*m 

Tjo.r derioi'ctj tlio conjugate complex , f(t) account c for 

the foroGO on the boundary. In tho prooont caae f(t) = 0. 

So far, the origin was talcen at the centre of L. How 
a mapping function T,^ y^(S) ^ ■i-c p ^ I? c'" is applied 
which shifts the origin to O’ , the centre of . 

'the transfonued equation of 1*0 KI=R • 

I'he bounda 3 ::y coaidition (47 ) is transformed to 


c|> Cl-) -+- 
{ 1 ^ 

'U> f (r-) 


-i- 


where 

c^, (.-) 

" 1 0 

■r: 

4^0 j 


4:'.Qt-') 

ts 




■= 

to 

^ lit 

and f - R-<- - 


on 

^0* 

functions analytic in the 

region )-£KR . 


(Ehese functions 4*1 found out 

using the methods discussed in chapter I. 

The expressions of sectioimlly holomorphic functions 
4>, (:$)=* and 4^, «)* 4^ may be written using 

(37) and (38) . I'or matrix they are 


ZfL c R- Cl -5) 

■d" 

CR-Cii)^ 


(Ai-cAO e, 

I R - C) •«) ^ 



g C( ^ ( 49 ) 

for li:-c| >Y'f 


(A( +t Ai) 8. 

Hi ll ■ ii m il iii n . m ■— 

W- C) 



axid 

^ Be, C:^R- __ 

' 2 CR-C|S) (|^„c,^)2- [-1-0)3 

[d - c,i;2- (R-e,l)’i 


( 50 ) 


[l5-c)>Y-' t 


where the following symbols have been used 


B= - 

(K +1) R 


+ iA^ =: 
1 2 


- Sg -f iCdj^4- 6g) 

(K + 1) R 


(61) 


h. - iD„ « - - eg- i(6^+6g)jr’'^ 

1 <5 ■■ '■■■ ' ^ ■ ■ n 

(K + 1) R® 

and 0 ^ *s o/R • 

From (49) and (60) and th® result given below in (62) 

p •+ P :z ^(ur 

fp V-V- a I _i , 

(68) 

jrtt?' 

Pv^^ - P,p W = * *' ■* '*'''“'11 . 


stresses for t^ie matrix may be calculated. 

It may be verified that the normal stress P^p and shearing 




stress vanish on the boundary of the circle HlsR , 

as they should* 

Uhe sectionally holomoi^hic functions <|)( 2 ) and are 


B (^R(z+c) + o^(z+c ) (A^-.iAg)o^(z-fo)^ 
2 R(E- c^(z+c)j (R-c^(z-t-c)J ® 

U^-iA^)(z+c f (D^-lBg) (z+o)^ 

R ^E^c^(2+c)^^ |BUc^(z+o)] ^ 


K+1 


Bo^[gR - Cj^(z-I»c) j 
2 {R - c^(z+c)J 


2 cR(l.-iA^) 

..,..„A ft. 

-[R- c^(2+c)J^ 


(63) 


2(A^~iAg)(z+o) { 2R »Cj^(z+o)} 
(R - 0^(2 + c)] ^ 


3(D^-.i3)g) R^(g+o ) ( i^>.lAg)Rz 

(r- c-(z+o)} r*^ 


f 
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(|) (a) « “ B ^R(z><-c) + Q^(2-»c)^ j. (i^-.iAg)c^(E:+c)^ 

2 -c^(z+ c)] {_^c^ (z + 


(A^-iAgXz+c)^ (Bl-iDg)(z+c)® 

R {,R-.c^(z+c)j^ {R -.c^(z+c)| ® 

(A^+iAg) R 

'■" ■ ' ■ ' '■■■ f 

z 


Bo^ {gR^ g cR(Aj^~ lAg) 

2 ^R- o^(z+c)5 {B -.c^(z+c)j^ 


(54) 


+ ^(y^yCz-i-c) -tgB-. c^(z-fo) } 5(l)^-.i33g)R^(z-<-o ) 

•(^R - Cj^(z + c) ^ ^ {R~c^(z+c)j^ 

_ (D^+ IBg) E® 33 

^ 3' ■** T * 

z 

Stresses in tRe polar coordinates (r, e) (z es r e^®) are 
determined by using tRe formulas 


^rr + ^09- 4 E® [<>'('=)] 

Pgg - ^rr ■*■ ^re ” ® 1>’'(a) + 't-'C®)] • W 
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Therefore , we gei; 




g In^+^ciz+o)--. 0^(z+c)^j 

2 


2 R ^R-o^ (z+c)| 


3(D^-lDg)R (z+c)^ 2(A^-il2)c(z+c) 

^ » 

^R-c^(z+c)} ^ ^R-c^(z+c)J® 

( A^~lAg)(z+c)^.(gR - 0 ^ (z4-o )} 

R |R ~ c.(z+o)}^ 


^{Ql + Gg -1(6^ - 

mm Iirm .-nrr niT- TIlTn . “ ■ — -f" II 

K + 1 

( 66 ) 

and - <^rr-^ ^ ^ He p|e,^eg^lCS,^6g )j 

K + 1 


+ dAg) R 


z 




Also 


(^00-^pr+2i ^re^l « 2 e 


2ie 


2Boz 


B o' 


jj nr g 

^R-c^(z*t“o)^ R(R«o^(2+o)|- 


6(B -.iI>g)R{^R( 2 +c) + o.(z+c)®} z 


+ 12 


(r - c^(z + o)} 


6 
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5(B^^iDg) R^{R H-gQ^Cz+c )^ 2c(A^~iAg)g {R+ec^(g>fo ) } 
'l^ 


{R-.c^(g+c)^^ «fR-c^(a+c)^^ 


6CA^-iAg)R(z+c) 4R^(A^-iAg) 4c'='(A^-.iAg) 

|^c^(z+c)j^ -^-c^(z+c)} -Jeu. c^(z+c)j- 


3 


(A^ - iA2)R 


..2 


(57) 


^re’m = +Pe21er(A^)E 


r 


• 2 


2(A^+iAg)R 3(l)^+iB2)R‘- 


z 


z* 


Z' 


Oontinuity of Bormal stross and, shearing stress P^ 

at tile equilibriiua Interface may be verified using (66) 
and (67). 

The displacement components in polar coordinates 
are given by the formula 

2^(u^+i Uq) « e“^® [k <[>( 2 )- z <)'(a) - • 


So that 
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{H- Oj^(z+c)] 

B z ^E^+Sc(z + c)~ 0 ? {"z +C5)^j 


2R - c. (¥ +0)} 


Bc^|3R -.c^(z+c)] K(A^-iAg)c^(z+c)' 

- ^ 2 

2 |_Er-c^(z+c) 5 -{R- c^(z4-c)J 


K(A^-iAg)(z+c)^ 2c(A^+iAg3g(z >fo) 


R|_R-c^(2+c)| {R- o^i z+c)|^ 


,+iA„)2(z +c)^i 


r(R-c^(z +c)}' 


•c^( 2 +c)| 2(A^+iAg)cR 


{R- c^(g+c)} 


2(iij^4-lAg)(z-<-o){2R - c^(z4o )| K(l)j^~lI)g)(z+o) ' 

|R « c^(z +o)}^ ^R- c^(z+c)j^ 


3(D^+iD2)R z(z+c)^ 3(I)^+iDg)R^(z+p) 

^R- c^(z + c)J^ ^R - c^(z + c)j 


>^K jj_e^ * ^2 " ^ 

K + a “ 

K + 1 r’® 






2;U.(u^+iUg)^ 




g5:j^ei+e2-i(6i~6g)j ^ 


^ (A^+iAg)iai /i-ie^^+eg+iCfi^-dg) j a 

5 

(A^-iA2)lz (D^-.i3)g)R® 

z2 r * ® 'z 

It may be noted tbat the relations (42) may he obtained from 
(58)* 

If the outer radius of the matrix R tends to infinity, 
tbie results for circular inclusion in an infinite medium 
are obtained((6)) . On putting c « 0, the results for 
concentric inclusion are obtained . IChe results of half 
plane may be obtained by taking R - c =* d where d is finite 
and R and c both tend to infinity ((6)). 

In the above analysis th© elastic constants for the 
inclusion and matrix were assumed to be the same. However* 
if the shear moduli for both the regions be the same but 
Poisson ratios be different, then the problem may he solved 
as follows : 

Ihe boiindary conditions in this case shall be 
^(t) + t |)* (t) + '4^ (t) « f (t) on (69) 


<)'^(t) + t ))''^(t) + =* <)"('*;)+ t <j)’lt)+f (t) on I. (60) 
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K - t - v|/(t) = Kj^r(-t) --t vfCt) 

+ 2jU-g(t) on It (6l) 

where K a 3 - 4V and K^«3 - 4 \j ' for the plane strain 
case and K = and generalised 

plane stress case, and v' the Poisson ratios 
of the inclusion and matrix respectively. 

Adding (60) and (6l), We get 

2M,g(t) , . 

(t)**' (t) - K‘ <b“(t) ^ on L (62) 

K+1 

Kf+l 

where K’ « 

Prom (f»0) 

y+(t) . f(t) = - L?^) - - 1 

on I. (63) 


Ecg^uation (62) may he reduced to the situation of the fonfl 


I^(t) - ^(t) » 


2 

K+1 


(64) 


when ij)(z) f® defined as 

(|)‘^(z) « 1 *^( 2 ) and K* ^"(a} « l*"(z) • 


(65) 
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Having found it)'‘’( 2 ) and ^(z) and may be 

determined from (62) and tb.e problem is solved. 

Stresses were calcu3.ated numerically and graplas 
of normal stress, shearing stress and hoop stress both 
for inclusion and matrix at the equilibrium boundary 
have been given in the appendix following this chapter. 
In each case 6 varies from 0 to n • Two cases & 

and -Gg « G have been considered. As a typical case 
the radii of inclusion and matrix are in the ratio 1*4. 
let the ratio of the distance c of the centre of the 
matrix to the radius of the inclusion be denoted by C , 
Graphs have been drawn for six values of 0 :0.0, 0.5, 
1.0, 1.5, 2.0 and 2.6. Comparison of concentric and 
eccentric case has been done whenever the difference was 
significant. 

It is interesting to note that when 6g => G 
and the inclusion is very near the outer boundary, the 
hoop stress in inclusion at the equilibrium boundary 
becomes positive for some values of 0 and is nearly 
zero when 0 « 2.7 . Por Gj^ ** -Gg « G , normal stress 
shows an interesting behaviour when the inclusion is 
sufficiently near the outer boundary. 




iPPSwDIl (DO 0ILAP(DBR II 
























oaaPTEE III 


ELLIPTIC lEOLUSIOH IH A CIECULAR EECIOI 

This chapter deals with the problem of misfitting 
concentric elliptic inclusion in a circular region. 

The problem investigated may be stated as follows t 

Let a circular ring of outer radius R have a 
concentric elliptic hole. The semi-major and semi-minor 
axes of this elliptic hole are a and b respectively 
(Figure 2 p.^8 ). This medium is called matrix. How 
in this matrix, an elastic solid of dimensions slightly 
larger than those of the hole but remaining within the 
limits of proportional elasticity are embedded. The 
elastic field developed due to misfit in size is evaluated. 

We use the terminology of chapter II. As stated 
there the outer boundary of the matrix is denoted by 
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and tlie inner boundary is denoted by 1. !Ehe centre of 
is talien as tbe origin and is denoted by 0, In tlie 
absence of tbe surrounding material, tbe inclusion 
misfits by a displacement (e^^x + 6^y , Cgy + bgX ) with, 
respect to the origin at 0. 

If the displacement components in Cartesian 
coordinates at the equilibrium boundary for inclusion and 
matrix are denoted by ) 

respectively, then 

- (e^x 4- 6^ y) a g^(t) , 

( 66 ) 

- (GgJ + bgX ) a gg(t) . 

It may be noted that ( uj , v:J ), the displacement at 
the boundary of inclusion is measured from its natural 
state ; the state it would have in the absence of the 
matrix ; while (u^^ , v^" ) is the displacement at 
the bo und ary of matrix measured from its natural state, 
which is the state when its inner boundary is an ellipse 
of semi-major and semi-minor axes a and b respectively* 

The equation of is IsJ =* R • ^or a point on 1, 
z is denoted by t* Substituting (66) in (32) and (36), 
it may be seen that 
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g(l;) ^6]^+ 


1 

i 


( 67 ) 


liCij) » (6j^+02^ ^ ^ 


+ (e^- Gg+iCd^+Sg)} t dct . 

2 dt 

Uow putting tlie values of g(t) and li(t) from (67) in (39) 
and (40) and evaluating the contour integrals, the 
expressions of f)^(z) and ^^(z) for inclusion and matrix 
come out to be 

h { \ iL- e 2 -i( 6 i~ 62 )}a - JL. 

V*iKz) » - (K+1) (a+h) 


(K+ll 


(68) 


jJl ^(e 3^+62) (a-b) 

“(i+T) ^ (a+b) (a+b)2 

2 2^^^^ 

2^ a b ^ 6^-62+1(61+62) ^{.2- (2 ""^ ) i » 


•- (a2«h^) 


1/2 

4 ^^ ab (61+ e2)^z-(z^-c^) ] 

^ (^2) 


(69) 


2 2^1/" 


2JU. ab(a^+b^) (ei-e2+l(6i+62)}{2 (z-(z -c) ) 

«2 

® ) 


2 ^2n2 


(K+l) (a^-b'^) 


(a 
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wliere c = a - "b , 

Tlie ‘boundary condition satisfied on the outer 
boundary is given by (4l). In the present case f(t)=0. 

Using (69), f^Ct) can be found out. Hiie value of 
fpC'^) is substituted in (4l) and then the functions 
t|)^(z) and which are analytic in the region 

|z[ 4^ R may be deteimined using the methods given in 
chapter I. Bie analytic functions ^^iz) and are 


t)o(z) 


Jd. 


2ab 


(K+l) (a^b^) 






2(^) I ^j. 1 

(a^-b^) 

(a^+b^) (ei-62-i(6^+62)}( 


and 


. 4 . / "N ^ 2ab i 




( 70 ) 


E' 


8 


(E^-0 zb®/" 


i mm 2(6^*!* ^2^“^ 


R' 




2„2 



- E^- 


.2( a^+ b^ ) . g e2-i(6-+ ^ 2 ^-- 


^2 2 
c z 


} ig + (a2+ 1 6^-62-1(53^+ 6g)} 



The sectionally holomorphlc fimctions (j)(z) and (z) may 
be found out using (37) and (38) and are given below. 


<)i(2) * - 

tl(£=) = t(z) + 


;JL 

K+1 

/jL 

K+1 


.[Gi-^eg- 1(61-63)} 
2 (ei‘^So) (a>.b) ^ 

(a+b) 




(VI) 


(63-62-1(63+62)} 


(a~b)^ 

(^)’2 ^2^! 


z , 


(tjz) « ^^iz) + 


2Ai 


(K+1) 




(72) 


2ab jU, 

14(2) = ^g(z) - 


2(63+63) (z- (2^-0^)^^ 


^ (63-e2+iC53+*2)U2(*-(^"-c^)^j- 2^3/2( 

(a -b^) (z-c) ^/di- 

stresses in Cartesian coordinate system are determined using 
the foimulas given in (9) . They are 
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4 ab 


(P + P - 
XX yy^l 


JJ.U 

(ej^-eg-i(5^+6g)}| — -^ 2 s /2 - ij 


2 (ei+ e 2 )R^|- 


,2 f __ R^ 

2 / -04 2 2 x 

2 CR -c z ) 


1/2 


1 c 

+ — + 

„2 4R' 

z 


4} 


Sis!i|!) E 2 ^e,-eg-i( 6 ,^ 62 )l 

(a ~b ) 

>2 


!• 


R‘ 


2 /T 34 2 2 \ 

z (R -c z ) 


1/2 


i + S-. 

^ 4 R^ 


1 


.(aW)^ei-eg-i(6^+5g)^-^-— 3/2 - ? 

(R ~c z ) 2R 


r :" 


+ -[S ^- 22 + 1 ( 6 ^+ 


R^ 


(E*_o® z 2)^/2 


■} 


2 (ej+ eg)a® j 


E“ 


z 


" -i Z4. 


2(H^-c®S2) 4E' 


H 2 (e,^g+l(S, 4 i 2 )^ . 


R"^ ./« + *■ + 


^^(R'^- c ^ zS ) 


1/2 


z ^ 4 R'^ 



S5 


• fctr i "{Sb*) > 


r> — “ O I 

2It ■' 


(73) 


= ^^xx-*-^5y^i ■" EFi{ 61+60+ Cei-eg)J 


(a-b) 


4a‘b 


^ (e.-0^+iC6.+ 6 p)l { 1- } 

*~*o -i- C > t OOi /O 

(a^-B") (z‘--o‘')V'' 


+ E+i , .£“ ^2 


. JL. {e,-e,-i(6i+5„)Kl- 


■S' 




} 


and 

, V M 48.1) 

(P — P, + 2 iP _)j= T^j.'i , 2”" 2V’ 
yy CT 3ty 1 n+l (a'^^b ) 




I o^z®+2®Vj) I _ a(e.+ So). 

/ ■rs'^C ^ -M ^ \ c ' / X iC 

(R - C Z ; ' 


I 


R'^(2R'^-3c^z^)2 

, '- pTi " W - T ,., 

(C *w \^ / 

3 yj.v ~C- 'J, ] ' 


'Z^z 


or 3# 


-2 


1 


2(a^+b^) 


^ R^(4c^s^-3R^) 

"^2? " ^ 

fe -e -i(6 )l I 2E!LI 

1 ® 1^2 ^^* 1 ^ 2-'/ 1 ^, 5,,.4 ^ 2 ^ 2 ^ b /2 


a^ va -b^'} 


Q 


^j-i. -“v..; 3 

■nS/ji 2 2 r...n4% 2 

R (4c a -oE ) c 

3'-^(E*~c^a^)V2 '*' 2z^ 


a 


4 


} 


a- 
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-(a +1) )(e3^-G2-i{6^+6g)]-( (g4_.^232)6/2 “ ^ 


B ^( 4 cV ^^) 1 . 


2 it 

K +1 


2(a— b ) (G^+6g)+^G^-Cg+i(6^+62)J (a— b)' 


(a+b) 


(a+b) 


+ (G^ - Gg - i(6i +^2^ 


(74) 


' yy X3C 3jy'ia ' yy xx. 3cy 1 — — 

K+l 


2(a~b) (e^+Gg) 
(a+b) 


+^6^~Gg+i(6jl^+62)^ (a— b) ^ 

(a+b)^ 


4. ;ii. ^ob 'L®1“^2'^^^^1'*^2^1 

(K+1) (a^-b®) 


L (22-c2)V2 


(a®+b®) 



1 


( z ®- o ®) 


1/2 



M- Sab 

( K + l ) ( a ^- b ®) 


( e ^+ Bg ) (1 


Z 


iM 



The displacemeirt coraponeuts in Cartesian coordinates are 
given by (lO). Ibr inclusion and the matrix, they are 
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2al) 


2 ja(u+iY)^ =“ fei *■ 


K^G^~€ 2 -i( 6 ^+ 6 g)]|(R^-c^ 3 ^) 


1/2 


- R^- 






4 r 2 


] 


i 

z 


-K(a^+b^) ^e^-e2~i(6^+6g)]|^^4^^2^2^l/2 " 


z 


(^1*^2 ■*‘^^^1‘*^2^}-^(j^4_q 2 ~ 2)3/2 ” ^ j 


+ SCe^+Bg) -( 


R"* 


- 2,^4 2-a ' z 

z (R -c^z 


R^ c^ 

Wo + ^ + ~o ^ Z 

1/2 4 r 2 


!l^5 ^ei-eg+i(6i4*2)|{- 


R^ 


(a®-D^) 




+ ^ +(aW)(ei-e 2 +i( 6 j^+ 6 a)] 


R^ 


^ (R^— 


21 *^ 


I z -(ejL-e2+i(6^+6g)j 


C^«4 .2 — 2%^/^ R® T 

^(R.c z ) ~ ^b 4^;,2-72)5/2 i* 


z 



+ 2(e+e)r— 5i— , -E^ sii®? i 


2 (a®+b®) ( £. _ „ rx Jj; ur . _ E* 


8)1/2 


}■ 






1 

z 


mk 

K +1 


'[^1^ e 2 -i( 6 ^~ 62 )j+ -.i^2.^ei-e2+i(5^+&g)J 


z 


JL 

K+1 


r 


{ei+eg+i(6i-6g)] + i±±i ^ei-e 2 -i( 6 i+ 62 )J 

(a+b) 


2 M- 

K +1 


(a-b) 

(a+b) 


(e^^+eg) 2 


— 1 ±^I^ 

K+1 (^b)^ 


Xe^-eg-iCs^+figy. 


z 


eg+ 1 ( 6 ^ + 6 g)] i 

K+1 


and 


( 76 ) 


2 jli (n+iv)jj^ a 2)^(u+iv)^+ M-K -(^®i'*'®2'^^^^l~^2^i+ 

K+1 ^ (a+b) 


*C® 1 ~ ®g ■** 6 g) ^ 


flllN 




p- ' 
K+1 

|e^+62+i(6 

1-62)} + 

(a-b) 

(a+b) 


.i(6,.&g)}J 

2 }l 

(a-b) 

(e-i+eo) 2 

+ M- 

(a-b)^ 

/ — So—i (6 

K+1 

(a+b) 

X M 

K+l 

(a+b)^ 

11 ^2 ' 






(i2-c2)V2 


i 


iL -ML- fp +R_W “ - V!J^l + ^ 2a-b(a +b ) 

“ 2 ^2 'z -c'^I 7“o — 2T^ * 

K+i a ~b*^ K+1 (a^-b'^)^ 


2 

{e^- e2-iC6j^+ Sg)} r 2(z - ( z 


It may be noted that (66) may be obtained from (75). 

The normal and shearing stresses on the outer, boundary 

L vanish which may be verified using (73) and (74)* In 
0 

order to verify "the continuity of noruial -and shearing 
stress across the equilibrium interface i.e. I», the 


transforoiation 

7. ^ C. COi-tv i 


( 76 ) 


> 
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is used ((4)) • The boundaiy of ellipse +y^/b^= 1 

is now defined by v/here ^ and c are given 

by the equations 

C CoS’L f ^ = O. 
and .f 51 -nfv^o = ^ 


If the normal, shearing and hoop stresses are denoted by 
^ respectively, then 

f ^ I 


and 


f 


•n-M 




.(TO) 

■S^ L T) S 


Using (73), (74) and (76), the continuity of normal stress 
and shearing stress P-^rj equilibrium boundary 

may be seen. The difference in the hoop stress at the 
equilibrium boundary is given by 

, (P . ._st 

i I (l<;+n(a>-U*-) • CoiZT]) 

( 79 ) 

Some particular cases of great interest may be 
derived from the results of this problem j for example 
if the semi-major and semi-minor axes of the ellipse are 



aiid tlio Gemi-niajor a:cis a of the ellipse has been 
the values 1,2,3. Por each value of a , the ratio a/b 
talces the values 2, 4 and 10, The ratio of 
shearing stress to /'•^/flC+l) has been denoted by T . 
nine {graphs for these cases are given in the Appendix 
fol3.ov/ing tills chapter. It is interesting to note that 
for slender ell.ipses, the lines of maximum shearing 
Gtref 5 s in matrix emonating from the boundary of the 
ellipse resemble to some extent with those of tliin 
rectangular inclusion in a circular region which are given 
in the Appendix to chapter IV. 
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CHAPTER IV 


RBOTAWGTJLAR INCLUSION IN A CIRCULAR REGION 

In tile analysis of tills chapter, the stress and 
displacement fields liave been evaluated for a rectangular 
inclusion in a circular region* Consider a circular ring 
whose outer boundary is a circle of radius E and inner 
boundary is a concentric rectangle whose sides are 2a and 
2b (Eigure 3 p.i^s ), The origin is taken at the centre 
of the rectangle which is also the centre of the outer 
circle. The outer boundary of the matrix is denoted by 
L^ and the inner boundary by L . 

The displacement components of the inplus ioti in 
Cartesian coordinates in the absence of the matrix are 
prescribed by (e^x, G 2 y) with respect to the origin at 
the centre of L^. 



At the equilitrixim interface 


and 


= - e^x = g^(t) 

T+ - V = - ®2y = 82^*^ ’ 


(81) 


where and uj , have the same meaning as 
before and they are measured from their natural states 
as stated previously . 

Prom (33), (36) and (8I), it may be seen that 


g(t) = - (e^+ 63) t - (e^- Gg) t , 

2 2 

C82) 

h(t) a (61+ G2) ^ + (61- Gg) t + (e^- 6g) t , 

2 2 dt 

Substituting g(t) and h(t) from (82) in (39) and (4o) 
respectively and evaluating the contour integrals, the 
expressions of (t)*(K) and v}^^^(z) for inclusion and 
matrix come out to be 


1)*i(z) » - /^(Gi+Gg):^ _ / ^(e^-Gg) 

(K+1) n i(K+l) 


z log 


/-(a-ib)^ 

z^-(a*ib)® 


+ a log 


(z-ib)^-a^ 

(z+ib)^-a2 


+ lb log 


(z-a)^+b^ 

(2+a)^+b^ 


J 





?5 



mmmrnniSTn i iii ii W m .f*. 

(K+1) 


g/^Ce^+Gg) 

Tii(K+l) 


z log 


z^-(a+ib)^ 


(2-ib)^-a^ 

+ a log 75—0 + i'b log 

(z+ib) -a* 


(2-a)^+ b^ - 
(z+a)^+ 


ACe^-ep) r (2+ib)^-a^ 

^ 2 1 a log + ib log 


Tci(K+l) 


(z~ib)^-.a^ 


(z-a)^+b^ 

(z+a)^+b^ 


and 


111*^2) 


TCi(K+l) 


z log 


z^-(a-ib)^ 

z^-(a+ib)^ 


+ a log 


Cz-ib)^-a^ 

(z+ib)^-.a^ 


+ lb log 


(z+a)^+ b^ 


S/lCe^+eg)^ z^~(ar.ib)^ 


(84) 
(z-ib)^-a^ 


'L*_(z) = £ z log —5 + a log ~~ :o ~ 0 

7 ti(K+l) L z^-(a+ib)^ (z+ib)^-a^ 


+ ib log 


(2-a)^+ b^ 
(z+a)^+ b^ 


^/l(G^“Gg ) 
It i(K+l) 


(&fib)^-a^ 

a log ^ 


+ ib log 


(z-a)^+ b^ 
(2+a)^+ b^ 


J 



She eciuation of is jz] = R. On the outer 
houndaiy , the boundary condition is given by (41). 
In the present case f(t) » 0. Using (84), f^Ct) may 
be found out. Ihis value of f^Ct) is substituted in 
(4i) and the functions and which are 

analytic in the region |z| 4^ R may be determined using 
the methods described in chapter I, IDhey are 




M (6^"G2^ r 


7ti( K+1 ) 


R^-(a+ib)^/ 




(R^+ib25)^-a^!2® 


+ ib log 


(R^-.az)^+ b^z^’ 


(E^+az)^+ b^z^“^ 


Sli'Ce^+Og) 


R^ log 

^ E^~(a-ib)^z^ 


+ a log 


(R^+ibz)^~ a^z^ 
(E^-ibz)^~ a^z^ 


it 1,. 1 


^eCz) 


^(Gi-Gg) 

1 xi(K+l) 


(R^-ibz)^- a^z^ 

a log ■ " ■ - T L. + 


(R^+ibz)^- a^z^ 


( 86 ) 




(R^-az)^+ E(a+ ib)^ R^a 

(R^4-az)^+ b^z^ (R^-(a+ib)^z^j 


2(a~ib)^ R^z 
^E'^-(a-.ib)^z^} 


2(a+ib) (a-.ib)^z 
{R'^~ (a-ib)^z'^i 


+ 


2(a~ib)(a+ib)^z ' 
{R^-(a+ib)^z^i - 


7ii(K+l) L 


R^ R^-(a+ib)^z^ 

? E^-(a-ib)2z2 


4;i,ab n 

""z J * 


The aectionally holomorphic functions (J)(z) and '4'(z) 
may now be found out using (37) and (38) an-d are given 
below. 


it)^(z) « <|)q(2) 


(K+l) 


/^(6^-8g) 

7ti(K+l) 


z log 


z^-(a-.ib)^ 

z^-.(a^ib)^ 


(z-.ib)^-a^ (z-.a)^+ b^n 

+ a log + ib log 

(z+ib) (z+a)*+b^ 


t (z) 


/i(ei-ep)z 

- to(2) + - + 




+ a log 


(K+l) 
(z-ib)^,,a^ 


(z+ib)'^-a' 


^ — -g + ib log 


7ci(K+l) 

(z-.a)^+ b 


(86) 

z^-(a-ib)^ 

Z lo g "h" ' p 

z^-.(a+ib)^ 


2.. -kS 


(z+a)\ b^ 


] 


+• 
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7ti(K+l) 


/^(Gi-eg^r (z+ils-)^ -a^ 

+ — — I a log - — — ^ g- + ib log 


Cz-a)^+ b^ 
(z+a)^+ b^ 


= (t>o(z) 


/l(e.-e2)r z^~(a~ib)^ 

iii(K+l) L z^-(a+ib)^ 


(z-ib)^-. vz-aj) +D -1 

+ a log + ib log , 

(z+ib) — a^ -t-h* J 


(z-a)^ +b^ 
(z+a)^ +b^ 


, ^ , 2JULCe.+ep;r 

%iz) = %Cz) + 2 log 


2JUL(e.+e2)r z^-(a-ib)^ 


(87) 


+ a log 


(z-ib)^-a^ (z-a)^+b^ 

( z+ib ) ( z+a) ^+b ^ 


/i(ei-eo) f (z+ib)^-a^ 

iti(K+l) L ® (z-ib)^-a' 


^ ^ +ib log 


(z-a)^+b^ 


(z+a3*+b 


2^-u2 * 


Stresses in Cartesian coordinate system are given by (9) 
Using (87) and (9), it may be seen that 


(P +P ).« - 
^ XX yy'i 


4ju(e-+ep) 2AMe--eg)r R'^.*z®(a+ib)^ 

mmrmrnmm-mmmmmmmmmmimmm ^ ILO^ mmm ™ mii Miiii n ii Wh i n i H iiiMii n in h i lu m 

(K+l) ■ , wi(K+l) L R^-(a-ib)V 


2R^(a^+b^) 

’{?-(allb)2/j*' (R^-Ca+ib)^z^}’’ tR^-(a-ib)^z^] 



2R^(a^+T3^) 




i;E^-(a+lb)2 ;2j * i°s 

R^-(a-.iTD)^ 2R^ 

R^-(a+ib)^ z^ (.R^-Ca+ib)^ z 


2R^ 

(R^-(a-i^^^ } 


2R^(a^+b^) 2R' 

(R^-(a+ib)^z^} (r:^. 


log 


^-(a+ib)^-j 

z^-(a-ib)^J 


4 j-L(6^+e2) " 

Tii (K+l) 



4iab R^ 

— K-~ + -p log 
R*^ z ^ 


R'^-(a^ib)^ z^ - 
R^-.(a+ib)^ z ^ - 


_ . . 4;i.(ei+eg) 

(1*^+ 


= ''xx-" ‘yy'i ’ 


and 




4^(6^— 6g) 
lti(K+l) 


(a+ib)^z z 


4;R^-.(a^.ib)^z'‘^} 


2 2 ; 


f, 

(a-.ib)^z z 2R'^(a+ib)®zi’ 

’*’ {R'^-(a-ib)^z^^J (#-.(ar»-ib)^z^|‘ 


'(a^+b^) 

.(a-ib)^z^l 

H^-(a^fib)^z^ 

R'^«(a-ib)^z^ 

( 88 ) 
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2 R ^( a - ilD)^zz 




2 R ^( a ^+ b ^ a - il))^zz 


gR ^( a ^+ b ^)( a + i - b)^zz 

+ ■■ ' " “ ■■ -nil r -, nr - T - ^ 

2 _ 2 ) 2 


z z 


{R^-(a+ib)‘^ z^J 


z ^-( a - ib )^ 


z z 


z ^-( a + ilD )^ 


2 R ^( a + ib )^ 2 R ^( a - il ))^ 

iR ^-( a + ib )^ z ^ i ^ \ R '^-( a - ib )'^ z '^]'''^ 


( a ^+ b ^)( a - ib )^^ RS ( a - ib )^ z^j 


(:?-( a - ib )^ z ^} 


2 _ 2)2 


(a^+b^) (a+ib)^ ( RV(a+ib)^z^] (a^+b^) 
+ -3 ,„.^2 . 2)2 * “2 


( R ^- ( a + ib )^ z ^|^ 


z - (a-ib)‘ 


(a^+b^) 1 4u(e.+e„)r 2R^z R^-(a+ib)^z^ 


z -(a+ib)' 


4 jU . ie ^+ Q ^) 
Hi (K+l) 


log 


z 


R ‘^-( a -. ib )^ z ^ 


2R^(afib)V 2R^(a-ib)^ z 

”** z [R^-(a+ib)^z^j z:jE^-(a-ib)^" z^j 

5R^ R^-(afib)^z^ 4 iab 

- — ^ log “x- — rrrr2 “* ”‘"2 — 

2i R ~.(a*-ib}“z z 

2R^Ca^ib)^ 2 R'^(a-ib)^ 

" B2{H^-(a*lb)®/} . '*' 2®(E^-(a-ib)Vi * 
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+ log 5~ 

z -Ca+ib}^ J 


2/1 (e^>- Gg) 
(K+1) 


( 89 ) 


P^J„ = (P -P +2iP ), 
' yy 2DC 3£y'm ' yj xx xy'^i 


SjJ- (Gj^-Sg) 

(K+1) 


The diaialacement components in Cartesian coordinates are 
given by (lO)* Por inclusion and matrix respectively, 
they are given by the expressions 


;|l(u+iv). 


__ 


M -(i’-K)(ej^-i-eg)z 

(K+l) 


>L(e,-e«) r p^-(a+ib)V 

Kz log -T s-p 

Hi (K+1) L R^-(a-ib)^z^ 


(R^+ibz)^-a^z^ (R^-az)Vl3^a^ 

+ aK log — s — —■"*" » a o + I'WK log •— r? — 

z^-Ca-ib,)® (z-ib)®-a® 

+ Ka log + aK log — 

z'^^-Ca+ib)"' (z+ib)^-a'=^ 


(z-a) W H^^i(a-ib)® f 

■*■ (zH-a)^rt2 " '•‘’® E4.(a+lb)® 5® 
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2R^!Z 

2R^(a^+b^)z; 


2R^ z 

R^-(a+ib)^ 3^ 

2R^(a^+'b^) 2 
R^-(a+il))^ 


+ 2 log 


? - (a+ib)^ 
? - (a-ib)^ 


- a log 


(R^+ibi)^- ,a^2^ 
(R^-ibI) ^-.a^ 


(R2-az)^+ b^ i® 
+ ib log 

(R^+ai)^+ b^ ^ 


+ 


2(a-ib)^R^ z 
R^-,(a-ib)^ 2^ 


2(a+ib)^ R^2 
El.(a+ib)^ 2^ 


2(a-ib)(afib)^ 2 
R^-(a+ib)^ 


2(a^+b^)(a-ib)^z 

R^-(a-ib)^ ? 


a log 


(z- ib)®- a® 

G +lb)^- 


+ ib log 


(2 -a)^+b^- 
(2 +a)^+b^- 


+ 


2 >l (e^+eg) I R^-(a+ib) V 

iti (K+l) L 2 ^°®R^-'( a-ib)^ 


(R^+ib2:)^r^a^2^ 

+aKlog . S ' : ' \S ’ ""'" 2 ~’ 2 - 
(R -ibz) -a 2 


O^-az)^ b^/ 
(R^+az)^+ b^z^ * 


+ 2(l«i:)ial) 




log 




+ 



log 


E^-.(a-ib)^ 

E^r(a+ib )^ 


4iab ^ 

"t:; — + z log 
z 


z^-Ca-ib)^ 


+ a log 


iz +ib)^^a^ 
(z~ib)^ -a^ 


ib log 


(z-a)^ 4- b^- 
(z+a) + b 


and 


(90) 


2 /X (u+iv)^ » 2 jx (u+Iy) j_ + 


( K+1 ) ^ (K+l) 


It may be verified from the expressions of stresses 
in (88) and (89) , that the normal stress and shearing stress 

vanish on L i.e. on Izl « E , which they should, The 

0 • • 

continuity of normal and shearing stresses may also be 
seen, provided the proper branches of logarithmic functions 
as given below are chosen* 

Let e^, Og, Sg 04 angles as shown in 

Figure A p* 48 , then for matrix 


I ( log 
I ( log 


z+a -lb N 
z-a- ib ' 


24-a>.ib \ 
z+a+ib ' 



» 


» 
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I ( log 
^ ( log 


gsMi]?.. ) 

z+a-hih ' 


z-a-t-lb \ 
z-a-ib ' 


and for inclusion 



I 


t 


I( log 


2 a - ib 
a - a - ib 


) = 0 . 


T/ n*_z + a-.ib 

5 " Va + li) ^ 


I( log 


2 - a + ib \ 
2' + a + il)' ' 


» 


I( log 


2 ~ a + ija \ 
a’"-' a -""i'b ' 



» 


where I moans the imaginary part of a complex quantity. It 
may be noted that the orientations of e^» e2» 65 64 

for inclusion are different from matrix. 

In the above analysis if the radius of the circle 
R ^ then the results for rectangular inclusion in an 
infinite medium are obtained ((9)). Ihirther, if a a b, 
then the results of square inclusion in a circular inclusion 
are obtained. Also if a « b and R - «> , the results of 
square inclusion in an infinite medium are obtained . 

Numerical work has been done for this problem and is 
reported in the form of graphs in Appendix following this 
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chapter. Lines of constant maximiam shearing stress have 
been drawn in these graphs for the cases when Ggs G 


and =5 “Gg = G 


IPor hoth the cases i.e. ej^MGg= G 


and e^s! -.Gg = 6 , square inclusions for a s l and a » 2 
and rectangular inclusions of sides a r= i ard ha .1 and 
a =s 1 and h ss ,5 have been considered. The radius of the 
outer circle is tajcen as 4 in all the cases mentioned 
above. Because of the symmetry only a quadrant of a 
circle is considered. It may be noted that the boundary 
conditions at the outer as well as inner boundaj^y were 


verified numerically also so as to have a check on the 
numerical computations. 


I 
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APPENDIX TO CHAPTER 


IV 

















CHAPTER V 


THREE-DIMEHSIOM AHE T^/O-^EIilBlISIOHAEj POIl-fT EORCE PROBLEMS 

In the previous chapters, inclusion problems 
when the matrix is finite, have Been solved with the 
help of the Hilbert theorem. However, when the matrix 
is infinite, energy methods ((lO)) and point force 
methods ((ll)) have also been applied. Energy methods 
((12, 13)) require considerable guessing of the 
equilibrium interface which is not possible except in 
some simple cases. Using point force technique, a 
systematic study for ellipsoidal inclusion was made by 
Eshelby ((l4)) and the two-dimensional complex variable 
formalism was first given by Jaswon and Bhargava ((15)) , 
The point force teclmique consists of the following 
hypothetical operations. 




Mrst, cut out tlie inclusion from the remaining 

material and allow it to undergo freely the non-elastic 

deformation. Next apply the surface tractions to the 

inclusion which bring it to its original size. Replace 

the stressed inclusion in the cavity from which, it was 

cut out. Rejoin the material across the cut. ipply a 

distribution of point forces on the inclusion boundary 

80 as to neutralise the effect of aforesaid applied 

surface tractions. This brings us to the required 

conf iguration . Thus, as it were, the effect of the 

deforming inclusion is to bring into play a layer of 

point forces along its boundary. Therefore! if the 

elastic fields due to a single force in an infinite 
medium are known» the corresponding fields for inclusion 

problem would be obtained by the integrated effects of 

the whole distribution of point forces generated along 

the equilibrium boundary. 

Subsequent chapters therefore deal with the point 
force, or concentrated force as it is some times called. 
It may however be remarked that point force finds 
applications not only in inclusion problems, but also 
in material testing, soil mechanics etc* 

In this chapter, the love’s solution ((l6)0 of the 
equations of equilibrium fox axi-symmetrio case in terms 
of two functions is described? secondly the displacement 




components for a point force in three-dimensional infinite 
modi-um and stresses due to a point force in two-dimensional 
region are given. For details of love's solution, one 
might see chapter (XI) of love's hook ((l6)). We giwe, 
first some important results of love's solution which would 
he needed in subsequent chapters for solving point force 
problems. Some theory is also given to make the thesis 
self-contained . 

Che equations of equilibrium for the axi-symmetric 
case (variation with respect to 6 is zero) in cylindrical 
coordinates (R, 9, Z) (x = R cos 0 , y = R sin0 , z » Z)are 


and 


^ ^RR 


P - 

. -^ER 

+ 

dR 

dZ 

•r '.Mil,, 

R 

5 ^ 

® ^ZZ 


8R 

9Z 

“1* 

R 



= 0 


=* 0 , 


i 

( 91 ) 


( 92 ) 


Ihe displacement component u^ in the 0 direction is 

identically zero in this case and the displacement components 

Ut^ and Ur, in the direction of R and Z are related to the 
K a 

strain components by the equations 


e 


RR 


ill 




e - ^ 

®00 “• T » 


a u„ 

“ O’ » 


®R2 “ 



®He “ ^ 


i 


Let 


®SZ = ° 

%z = - 


8^ 

dHSZ 


(9S) 


(94) 


where ijp is a function of R and 2 only. 
It may be seen from (92) that 


? re ra * 

22 g 


1 a j) 

R aR 


(95) 


no function of R is added, because any such function 
can be included in ^ . 

Prom the e3jp3:^essions of e-^ and e^g gi'ven in (93), 
it is clear that 


®HR 




(E 


'69 


( 96 ) 


With the help of stress-strain relations 

^EE + 2A*-®HE ’ 



4 


p 

■^ee 

A A +2 ^ 

®90 ’ 

^zz 

« Aa + 2 

®ZZ ’ 

^E9 

* 


^EZ 


» 

*P 

•^ez 

= 2^eg2 

f 

4 = 

^®RR ^ee"^ 

®zz) 


(9V) 


stresses as 

wlienoe y 

(I+V ) (Pjgj- Pgg) = s 

let a new function P(E,Z) be introdnoed wbioh is 
given by the equation 


(98) 


then 


p M + I* • 

'‘=^RR 5 2^ 

the equation (9l) can be written as 


(99) 


-||,+ C:Pge-v>^BB“ ° 

It may be seen from (lOO) that 


( 100 ) 


Pqq = 5 


(lOl) 



no arbitrary function of Z need be added, for any sucb 
function may be included in |) , ^ ^ denotes 


6® 1 _S_ 

a ® an 


d 

d z^ 


the subjects of operation being independent of 0. 

All the non zero stress-components have now been 
expressed in terms of the two functions (j) and f .. from 

the fact that the sum of the principal stresses is a 
harmonic function , it can be derived that J) satisfies 

the biharmonic equation 

R r dR'*'az^^ 

( 102 ) 

The functions ^ and f are not independent of each 
other and the relation between them may be derived as . 


follows : 

Ihe eauation Ujj /a = can be written as 


% = I { ^ee ■ ’ 

where B is the loung's modulns of elasticity. 
Also from (96), (99) and (lOl) 


can iDe written 


and then the equation = 2/ie- 


^R2 


a u 


2 


B 


may he written as 

auf 


a 


1 


dR 


(1-^) + 



.i. 

(i+o) ap 

"D 

(103) 

6EdZ 


B az 

®ZZ 

k 

^zz~ ^^EE'*'^ee 

>1 

^ + 

1 ^ 
1 "i! 

p - ^ t) ] . 

(104) 

and (104) 

are compatible 

if 

’2 * ^ 

1 

az^ 

» 

(106) 


If a new function SI is introduced hy the equation 

9ll 5 

EP = + , (106) 

dE as 

then from (105), it may he seen that 

— » =i (I-aJ )y ^{j) , (lOV) 

where , as before, no arbitrary function of 2 need be 
added . 

Ihe stress components are thus expressible in 
terms of the functions ^ and SI. as 



a? 


a^)) 

1 

a(b 

^ER - j, 20 ■*■ 

- ( 

E 

— + , 

aE 

Pee =^Jv^{)- 

I ^ 

dij) 

5E 

d^(j) 

E-pa 

anaz 

» 


d^(t) 

1 

d(t) 

*D ^ ■■■ - — .« 

^ 


8E 

^E0 ^ ’ ^ez " 

Oj 


dE 


) , 


(108) 


When is introduced in (103) and (104), they may be 
written aa 


dU2 (1+'^ ) 


E 


5 dSL 


aj) 

»r ^ 


and 


au. 


z 


dZ 

Hence 


(l+'t) ) 9 , ^ ) . 

^ C Fz ^ 


B 


(1+i) ) , ^jx ^ ^ . 

^ - - _ — * aR ^ ’ 


Ug « 


c It - H-) • 

B 


(109) 


!I!he 


value of A can be found out in two different 


ways*' ' 

firstly A » ©HE •*’ ®0e ®ZZ 



A3.B0 


wliere 


A 

e 


9U.'p 011 r/ 

dF^" + -^ + ^ (110) 

n v^4 - ' 

^e:. 0ii:^ui±2^ ^^4, ^,^3 

EE 

is tiie sum of three principal stresses. 


® ” ^RR ^00 ■*' ^2Z * 

It follows from (llO) and (ill) that besides (10?),A7. 
also satisfies the equation = 0. 


Thus 

and 


= o 

3^ — ^ \7^4 

5x1- 


( 112 ) 


This completes the Love's solution of the equations 
of equilibrium for the axisymmetric case in terms of the 

functions J) and -H- . 

We now give the displacement components 
due to a point force (X^, 1^, Zj acting at a 

point (x^,y^,2i^) in a three-dimensional infinite medium. 




(Aty) ^ 




+ 


SIT 


•jLX.(.X'ty) 




( 113 ) 
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2 (x-h) ^ (y-k) 
{(x-.h)^+(y~k) (Cx~h) ^+(y-k) 


(l-Et)) (y-k) 

:s~: 


2)2 



471(1- V ) 


(114) 

(l-2^^) (x-h) 2(x-li)(y-lc)^ ~ 


The above stresses can be derived from the complex 
potentials given in Green and Zerna ((l)) or from 
Muskhelishvili ((2)). 



lOl 


CIiiP!DER VI 

ELASIIO CONE UNDER M M1£L EOROE 

Ihe solution to the problem of a point force 
acting along the axis of a cone at its vertex is knovm 
((l6)). In the present problem, the axial force is 
assumed to act at any point on the axis of the cone. 

This problem is physically very interesting. The solution 
is obtained by the following operation : 

Imagine a point force (0,0, Z^) acting at a point 
(0»0» d} in an infinite elastic medium* This states that 
a point force of magnitude Z^ acts along the Z-axis at 
a point d from the origin. , The expressions for 
displaceiaents are known and are given by blie 

equation (118 )i from where the stresses everywhere in 
the elastic medium can be found out . Ihls we call the 
first stress system. How suppose that there is a cone 
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vori;c;.: tlio ori^jin and the axio v/ith 

tho . Cti i:]u; ('-iirvod uurfrioo of the coiig, v/e can 


find tlio r,t:L‘Oi.i;i 


v/e can find the 


oji the firct utroc.ij 

noriiuil and ohcarini; etro 


ey iJ t eKi . Parti cnlarly 
rjfjes on the curved 


> 


Giirfn,c !0 of tlio cone. 

V/o now eolve another elect ic prohlem by applyinii' 
a oucoiiti, c^o.tom of force:.! on the surface of the cone. This 
s.'Ctvnd .Myo.tom of force:: are obtained from the normal and 
uhnarliii: sti'os.jo;: opiuosito in ;:i^’n to those obtained from the 
fir.'it stress system on the aurfo.ee of the cone. Hie tv/o 
stress systems in the cone are then superposed. Ihe solution 
to the problem when the curved surface is free from external 
forces and the cone is su.hjected to an accial point force, 
is then obtained, 

Gonsidor a semi-infinite elastic cone v/ith verter 
at the origin and the a-ac:is coinciding with the azis of the 
cone. She somi-vcrtical angle of the cone is 0^. let a 
force of magnitude act at a point (O, 0, d) along the 

a.>:is (figure 1 p.139)* 

Ihe displacement components 

0 arte: 3 ian coordinates for the first system of stresses 
whan the medium is supposed to be infinite are obtained 

from (llS) by putting and to be eero sad so also 

They are 
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t3wjlL(A+i^) 

(A-+/A.) y(s~d) 

ti.tttmTiiti|.Trl.irr-..i."4ii4m-i mi-Ti ■'i.ri -ni. m ■ 4 ...U ^ 

ulC j-*- ( ) Tg 

(A+/J~) z^ (A-t3^)2^ 1 

Q%lj.{AdZjx) rg ( A+i^) 


(115) 


wlioro r,, - x + y +• (z-d) . 

Af 

We ULio aphorical polar coordinates (r,e,«f ) so that 
X Rs r Bine coscj) , y = r sine sincj) , z = r cose (O ^ r < <», 
O-^^0^7T,O^'?< 2%) , It may be noted that 0 in 
cylindrical coordinates is different from e in spherical 
polar coordinates, Ihe displacement components in (115) 
may be written in spherical polar coordinates as 


w 

r 


(A+M-)2n (rp«d}(r-dp) _ ( A+l/t p 

8itjjL ( A+'a/^ ) (r^+d^-2di^)‘^/'=^ 8ii^ ( A+2ya ) (r +d -Sdrp) 


:i/2 ' 



(A4ja)Z^ dq.(rp-d) 

il» » V '» ,4» | l f M l >»tr ! rWH' *TTr.l IHl 1 **"* lO 

But jx (A'+A/i.) ( r +d -2drp) 


(A-tS/x )Z^ q. ^ 

(r^+d^-2drp)^'^^ 


( 116 ) 


tt 0 , 

where p ~ cose and q = sine * 

Ihe corresponding stresses are 
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z 


rp 




)X z 


0 


d(2p"^-l) 


4u(A + ^/i) (r'^+a -2d2p)* 


3 CA+A<-^ Zq (ip-d)(r~dp)^ 

4:% (r^+d^- 2 drp ) 


jX Z^ ^(r-2dp) 
r0^1, 47t(A+^yji) (r^+d^-Edrp)^^^ 


■ 5 (A- f At) Z^ 
411 ( A+ 5 ^) 

4 it (_A+ S/^) 

3 ( Zq 
4it (A+i/^) 


^dpq^r^- d^qL(l+p^)r +■ d^pg^ 

( Ad^-Ed 3 ^)^^ 

(pr + d(Eg^~l)^ 
(r^+d^- 2 drp)^^ 

d^g^(d-ip) 

__ ^^2 , 

(r^+d - 2d3^) 




/A ^0 

4it 


(rp-d ) 

(r^+a'^-ESip)^/® 


(liv) 


The subscript 1 is used to the above stresses 
to ihdioate that the stresses are due to the first 
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stress system. 

We now apxDly a second system of forces on 
the surface of the cone which are obtained from the 
normal stress (^00 )g shearing stresses aaid 

(^0 )g wliere , 

= -<^6e>l ’ (^re>g = a^d , 

on the surface of the cone. 

Hence 



SS — 

e « 0Q 


■ A \ ( Ppr -j-aCEgg - 1) 1 

_ 4ii CA+ ayuL) ( r^+d^-2drp ^ 


+ 


3(A+yU-) 
4Tr(A+ a/x.) 


d\o (d-iPp) 

(rSd®- 2drPg)®/®_ 



4ll (A-fSyU) 


qp(r~2dpQ) 

(r^+d^-2da^Q)^/^ 


5(AtA) Zq (I+Pq^ •«-a^PoqLo ;-| 

4n(A+a/x) (r^+ d^- 2drpQ)^/^^ J’ 


/p -N -s 0 i where p = cos 0^ and q_^=: sin 0^ 

' Sf ^ 


C1I8) 


f 

!lJhe mathematics tequired for this auxiliary problem 
is not trivial and is interesting. Some details of 
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calculation are given in the next chapter from which the 
solution has "been "built up. At this stage it suffices 
to mention that the second system of stresses 


come out to be 




^ 0 ° B^^cos mo; - ^ 


__ 

47C o E (l+lo^) cosh %a. 


(P, 


e 0'2 


471 ° B (l+ 4 a^)cosh Tta 






f 


00 EggCos met - E^^sin ma 


66 ^ 


471 0 B^^Cl+da^) cosh Tia 


da» 




(At^> 


00 Eyi^cos ma - ^ 88 ®^^ 

f r — ^cc» 


471 0 cosh Tca 


(P 



» 



and (^0(^)2 ° ’ 


( 119 ) 


where m^ «? log Cd/r)j ex.i® the variable of integration 

and the q,uantities. B^^ , ^2* ®33» %4» ®65’ ^66* 

and Bqq are given below , 




v/licrG ill the aho'VG e-rprocaion tiii'i the ercpreoiiione belov;, 

P (eofj o) ifj the TjOi^:ondro function of first hind and 

i ! 
r?* "t 

p’ (cooG) in its derivative v/ith roepect to tlie 

- ir + ia 

arijumont in the 'braohet; 


^11 “ ^^2’ 

= Hp + 2a Ih , 

®33 " il5 - 2a 

^44 " + 2a Hgj 

^55 “■ ^5 ^6’ 

^66 ” ’ 

Eijf^ = *“ 2a f 

^88 ~ ^ * 

where 

t\ « - {q.^(l"2V ) P‘ (p) + /( I + a^) :2 (p) } ^1 + 

- ^ +ia - J + ia 
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+ a f 21 ® P' Cp) + p( t + a®) P (p) 


S 


+ioc 


^ + la 


E, 


+ P’ (p) - 2p ( l+a^) P Cp) j Eg 


1 . 

^ +ia 


+ ia 


2 aq. P’ (P) . E 4 , 

- ^ + ia 


|^q.^(l-2v)) P’ (p) + ^ ( 1“ + a^) ^ (p) j Eg 


- ^ + ia 


^ + ia 


- a (2q^ P’ (p) + pC 1+ a^) P (p) 

^ -w+ia - 4+ia-^ 


l^P' (p) - Sp(|+a®)P (p) |B4 

1 . .. .. i 1 ■; - J 


- ^ -f-ia 


g + la 


+ 2 a q.^ E’ (p) - Eg , 

- i + ia 


p^( -2.5 + 2^ ) P' (p) • 

- l+ia 


!• + a^) P ^Cp) 


E. 


+ia 


afp®P' (p)+p{|>a®)I’ Cp) jBg 
i +ia - ^ + ia 

2 . o_2. 


2 

i . ^v 


+ (2p C P iCp) +(ciW) P\ (p) • B3 + 
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+ S a P‘ (p) • , 

- I + ia 


fp'" (-S*5 + 2^J ) P* (p) - % ( I + a^)^’ (P) I Eg 

1 L,- . ‘ - ^ +ia 


- § +ia 


+ a f E* (p) + pC I + a^) E (p) 

i- ~ +ia - ^ +ia -' 


E. 


+ f8p( |- + oe®) E (p) + (Sp®+ a®) E* (p) \ 

L - I + la _ I +ia J 

- 2 a q^ P' (p) . Eg ^ 

- ^ + ia 

« f(-^\ 2-2-^ +\)q^) P’ (p) + p(l-2^))( J +a^)E (p)}e^ 
^ i +ia - 2+ ia 

+ a (p^+ 2-0 q^) P’ (p) ' Eg - 2 P’ (p) * Eg , 

f-« ^ •f'iLoc — • ^ + i(x 

„ j'( E.\ g _ go+^),l^) E' (p)+p(l-^^>)(| +o;®)P (p) ISg 

I -| + i“ -l+iaj 


a(p^ + 2"^ q^) P* (p) • E^ ~ 2 P' (p) • \ , 


- I* +ia 


V& 


- ^ +ia 

1 


/2(l-V )+ J +a^]p P’ (p) + ( ^ +a ) ^ (p) 

- l+ia " . 


+ia 


\ - 
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I + a®) 5 (p) lEg 

“ ^ + i-a 


+ 2p a E' (p) ■ \ 

•It ^ 

- 'K- +ia 


% 


(2(1-1)) + J + /}p (p) +( |- + a‘) E 1 


1 . 2. 


- t+ia 


- 2 


p P‘ (p) + 2( i + a^) I* (P) 

A **1 

■i • 4 


« # + ia 


- ^ +ia 


E. 



(P) ] Bp 

+ la 



Ill 




+ia 


+ia 


■Po[4(l-%) )®+l -D +( 4 + a®)g®{ I - 2^; )}P' (p ). 

~ ^ia 


.P’ ^ + 3 ( 1 ~’^)Pq( I + a^)P* ^ (p ) 

- ^ +ia ° ° ^ - i + ia ° 


1 - U +“^>( I +(1-2"’ )(1+ k So )] • 


P’ 


(-Pq) :e 1 ,, (p„kP<,( 1+ a")®( I - S'!) 


- ^ +ia 


+!■« 


. P - (p_) P - /'•«"> 

- I +ia ° - I +ia » 


a 


S(l--0 )Po ( I + o^)P' 1 1 


+ ( I P^(1-2V )]P- ^^^(-Po)J_ 1 


Po( i +*^) 


- J+la 


(Pq) ^ 1 * 1 

•i» w + l-Ct 




We nov/ superpose the two stress systems given 
in (ll7) and (119) and obtain the stress field in the 
cone . It is 

, 

^oe ~ ^^00^1 ^^00^2 ; 

d ^ (l20) 

^re ^^re^ ^^re^2 ; 

p^^ * 0 and Pq^ =» 0 . 
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GHAPM VII 

BLitSIIO COME UEDER M mX1& EORCB (COETIiroED) 


In. tMa chapter the mathematics through which the 

strooKeo due to the second system of forces are obtained 

as discussed in Chapter VI, is described. 

Ihe equations of equilibrium are given by (91) and 
(92). lEheir solution depends upon tvm functions J)(R,Z) 

and »n.(H,2!) as stated in chapter V, Also as stated 1here, 

(|) and M satisfy the equations 


V-il 



i -i. 

'R dR 



= 0 


( 121 ) 


and 



(1-^3) V^<t) 


( 122 ) 


and the corraspo^JidinS stresses are 


given by (108). With 
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tlio htjlp of tiio 8trcB-iQB given in (l08) in cylindrical 
polf:'X Goordimitca , we can find out the stresses 

^ ^ spherical polar coordinates. 

I'Ms may be done by first transforming the right hand 
sides of (l08) into spherical polar coordinates and then 
transforming the stresses in cylindrical polar coordinates 
to spherical polar coordinates with the help of the 
trfinsformation formulae given in ((16)). 

The stresses in spherical polar coordinates are 


(]? 5 


8^ 2d 2p d 

( ^ ^ ^ ^ ^ 
r 0p r dr 


&P 




S. h y , 

r '“dr 


E ^ )_a , 
r dp ' * 


^^soh ” r dl ■ ^ 
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® ) <) + I (r-® 
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r 'dr 


, 

r dp ' * 


« a I ( i "■ K 


V j. - ^2 dp 


p d P d 
} if + 5 ( 6? - r i§ 


5-^}, 






+ i -2. 

r 5r 




) <t) 





drdp 


4 gS (W) g 
+ __ _ 

r dp 



j_a 



f(r,p) 


Auuiuiiinii; thtit 


V +ioo 

T . 3?"'® f(a,p) ds. 

Y -lOo 
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(lE?) 


8^ (b.jn.) and r® 8 , 

KV, na r _8 ^ 

OP Qp 

fU3 p “*■ j ^116 following relafioiiB may be 
ea.'iily obtained. 


J **Ay(it),'^)dr « s(a+s)(t) ,-^i) , 

0 ap'" 


J r® (^,-^) dp a - a( <j) ,-£^) , 

u>. M Jb 


dp « ( , 


(128) 


/* r®*"^ 


• • _ * « 


(i|),di) dr » ((j) , -H. ), 


wtor® daeh. denotee differentiation with respect to p 
and bar denotes the Mellin transform of the function. 
Multiplying equations (123), (124) and (185) by r®"^^ 
aM integrating with respect to r from 0 to «• and 
maJciag use of (128), it may be seen that , 


g.^ f*’ - 2(p <!>’ +s 1) - <i^(P'^'+ 





a a'"' ^ P - p^(pj-L» + s ja ) , 

Wu ' t 

- P () - JI' ’ -p ( 2q^s + 

(4-3p^))ii’ + s((s+l) p^- (2-p^)}ill , 


-Cl|^(l+s) + pCpjfi* + sJ^)}, (129) 

aluo ci^ jtTJ * - 2pix* + s(s~l)J^ = 0, (l30) 

- Pq.^(3+^Q)-nl‘ + s |p^(s+2)-l}j^ 

« (l-.*0){q.^ - 2p ^’ + 8(s-l) . (131) 

Equation (130) is a Legendre's equation whose 
general solution is 

^ . I P^_,(p) + B . (1S2) 


where A and B are arbitrary functions of s, is the 

legendr© function of the first kind of order s-1 and 
, is th® Legendre function of the second kind of order 

s-l ♦ 

Subetituting (13B) in (l3l), we get 
( 1 _V) If" - Sp f'+o(o-l) H = ® ((l+ 2 B)P®-sp^_iP>+»^lP )5 


pq^d+Bs) («;.!(?) + Bq-_^(p)). 


(133) 



iiojxitioii of (l33) is 


^ “ (r-vi “ Vi'p) - I 

+ g P1^(^g_l(p) + . (134) 

v/hoi‘(; 0 and JD arc arbitrary constants, Legendre function 
of the second kind . (p) , together with its derivatives 
with rer.i'oct to p, have singularity at p-1 . Since the 
sfcrot,u-eu shuuid remain finite (except at the point (0,0, d) 
on the axla)» the equations (l32) and (l34) cannot contain 
Q .(p) unfl Q„ ,(p) . So that B « o and I) * o and the 
BiiuirtioiiB (13K) imS (134) arc roduced to 


ji n APg_^(p) (135) 

and 

$ . J- fo P, ,(p) . ^ A bp" P 3 _i(p) 4 I W P;_i(P)].(l36) 

U«-d) L ' 

On oubstituting (l35) and (136) in (129), it may be 


aeon fchnt 




gO(l+D) |'(P) + ®P ^3^P^)J ’ 

^:^p(2"'^ -2+s) P^(p) +8(1-3 )^gCp)] 

^ + 20 |^(sq^- P^) I’s(p) 8 p(1+b) ’ 




(S-Sv « sp^- 2\) aq,^) P*(p) 

B 



3p(l+s) ) Pg(p)^ + 20 P^(p) 


n 




Cl 


{(2~2^J-D«c;^) p P’(p)-. s(l+s)P»(p)J 


+ 2G(l+a) (p P»(p) „ a p^(p)J 


(137) 


I'hu oaxiiituniti A imd G are determined from the two 
boundary cHuitUtionK given in (ll8) of chapter VI. We 
multiply (ilO) by and Integrate with respect to r 

front 0 to « * 5?he infinite integrals on the right hand 
aide of (llfO are evaluated by using the result given 
belov/ ((18)) 


m 


f 


0 


^8-1 


dy 

(l*t2r coso+r 




S# 2^“ ^ (sine) ^ ^ +^)B(s,2^-s}, 

. P S " ^(0030) (138) 

a p g 


0 < Re (a) t < 0 < K . 

fher® am other reoults which are also used. These 
aro given in th® Appendijc following this chapter. 

It is Buffloiont to mention here that talcing 
-2 < EiCs) < 0 la the complex s plane, the integrals are 
all fomd to h® valid and the boundary conditions in (118) 



are t. to 


(r "* ^Uo) {-/^+ ('^•f)(s~l)] (i^P^(-p^) 

0 1 0 {i 


-|ip(l+s) Eg(-p)l , (139) 


”» 8* L 




whom 


,(140) 


II( n ) 


n d 


s 


4n:(A4yA,) sin ns 


On imtiinr’: tho vcCLues of (^^(^ 00 ^gl 

Crain (I3y) mici (l40) in second and fourth, equations of 
(i;§7) rind no 1 vine: thorn, we get the values of the 
connfcantti A and 0 us 


A A 


l%(s) ^|(Pot g) 




fillli 0 




4^(Po»®^ 

E;( s) 43 ^^(Pq,s) 

AjCPq, s) 


(141) 


(142) 




Hj^(a) » 2{1-'J ) lisi . 

>) - +^a+a)p„P,(p)Jx 

fW 'o +(1+8)^-/' aA^)*) «j 
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/^PoCl+a) , ( 144 ) 


Aj.(Pt,,s) = [{Ah.J!;i)P^P^(_P^)+ (i+g) Pg(-P„)J. 

+ [f-/‘+ (At/i)(3-l)j q® Pa(-Po5-/^Po(l+P^Ps<-Po)] • 

.[(2(1-1))- 3(l+a)} p^ P;(Po)- s(l+s) PgCp^)] , 

( 145) 

« (s(l«v) - ^^0(1+3)] Pq ^s^Pq) ^s^Pq^ 


- Sd+a) q® p;(p^) P3(Pg)- s®(l+a)®p^p^(p^)p^(p^) . 

(146) 

Substitutiixg tbe values of A and C from (i4i) imd 
Cl4l0 lu (l37) and inverting with the help of the 
invoreion fomula (127), the stresses in complex Integral 
form are given by 


•/ +1 


00 




rr-'B Bui 


/ 


H(s) r 


rS-2 ^ 


V*!^ .5 -d 3 (p ,s) L 


A(pQ,s)[q''(2v; -2-.2s)P;(p) 


+ p(l+8)(2+i)d Pg(p)]+ 24^(p^,B)(l+s){a^Pj(p)+spP;(p)jJds 





eo'B “ an 1 


r-f 


H(b) 


V-ioo sa^(p ,s) 


^l(Po»s)^p^(s~2+2^J)P^(p) 


+Pa(l~s^) P^(p)| - 2\(p^,s)|s(l+s)p J'gCp) 


+ P^) P*(p)} 


ds 


(P ) 

' ^yW ' 


I 


V 


+10O 


g.(g) 3^’ 


-s-2 


Slsi J A C A f ^ 

TT V -l«o S4j(p^,s) 

- 3(i+s) (1-a^) ) p P^(p3} _ a^^(p^,s) P'(p) 

(3?^) 


a, (P^,s)|(a-SD _sp®-ai)sq®)p^(p) 

ds , 


re^s 


Siti Y-i®® 


s-2 


^‘^3CPn»®) 


^lCPf,»s) 


ICs-a-o -s - s^)p p»(p) - s(i+s) i’g(p)] 

• 24 (p^,s) (l+s)(p P'(p) - S P fp) )1 ds , 


(147) 


^^0*P^E ** ^ *“ ® * 


She solution obtained in complex integral form (147) 
i@ not affected by the particular choice of V , as long as 
l@(e) ( saV) remains within the limits -S < Ee(s) < 0. 

Ihe choice of V « - l/2 , which in turn Implies 




+ ia 


(l4S) 


B a — 


- «» < Ct < 
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ii.5 su^::,';:oo*t;ed 'by 
^^(p) together 
B « - I + ia, 


the fact that the hegendre function 
with its derivative is real valued for 
as may be inferred from the result 


E 


(cose) 

+ia 


2 

% 


cosh n a f 
o 


cos ax dx 


(Ecosh x+Ecose)^'^^ 


-71 < 0 < TC . (149) 

In order to establish the validity of our choice 
of V , we have yet to show that the integrands in (14?), 
regarded as functions of s, are free from the singularities 
if (l48) holds. The singularities of the integrand may 
arise from the polos of H^(s) and zeros of A-i, (Pq»s) . 

The polos of K^(s) are, in fact, the zeros of s sin n s. 

It may be noted that s sin sn does not vanish for s 
given by (l48) . 

It may bo seen that Ai,(Pq, s) is also free from 
zeros by considering the function h^Cp^, a), instead of 

s). ' Differentiating ^>(Pq, a ) with respect to 0^ 


8 43 (Pq» a) 
8 00 


a D ^ 

% - ^ +ia 


(Pq) [( 1 +i>)( I + (p^) 




Cp ) 

+ia 


, 0 < 8q< % . 


( 150 ) 
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llBindJ tto.e result 


•P “(COSO) = (-1)^1- 

P d(oos 0 )“ 

aiid in (l 49 ) t it may "be concluded that 


^ 1 ,, > 0 
~ +ia 


Oonsoquontly, (l 50 ) yields 


(for 111 = 1,5 -co<o^<«)) 

(I 5 l) 


a 43 (pQ>a ) ^ ^ 
90 . 


0 < 0 ^ < 71 
0 


Also as 0 Q 0 


a) «- I (l+V )( J » i*e* < 0 . 

Henco, it may IdQ concluded that a) ^0 for 0 < 9 ^^^ n 

and -.« < a < *” • ® ' I 

Beparato the integrals in real and imaginary parts. The 
real and imaginary parts are found to he even and odd 
functions of a, respectiTely. Oonseciuently, the contribution, 

of the imaginary part is zero and the stresses (Pp 3 .)gf 

N ^ (P ) . (P.jo ^ obtained 

as given by ©equation ( 119 ). 



1^5 


It, may be noted that while solving this 
problem of an axial force in the cone, another auxiliary 
problem when the tractions are prescribed on the surface 
of the cone is also solved. Hence all those problems in 
which tractions are prescribed on the surface of the cone 
can be tachled with the help of the method given for this 
problem after suitable modifications . 

In the next chapter, the convergence of the 
integrals in (119) is discussed. lo have a feeling for 
the stress distribution, graphs and tables are given for 
some particular cases. 



iffPEEDIX TO GECiPTBR VII 


The following results for our use can be 
deduced as particular cases of the result given in 
(138). 


s+3 j 

r r d r 


% d 


s-1 


(r^+d^--2drp)^/^ 


— ^ j- (2+sa^)Pi(~p) 

Sg'^simts ^ 

+Ps(l+S)pg(-P)j 


.c Rt -c 


oo ,8+2 _ 

5/2 


/ / ■ / ■T— .5/2 = e;(-p)-.(i+s)p^(-p)j , 


0 (r'^+d 3q sinus 




„s+l 


u d^**^ f /~2/„ .1 \ ov,2^■D^ 


/ f "p ^ ■ ^ —T S/£ = { (1 (-P) 

. Q A .. 3 q, sinus 


0 (r^+d^-2drp) 


dr 


+P s(l+s} PgC-p) j » 

-z ^ R< C..S) ■< I 


o (r'^+d -2drp) 


_jLi_ f_ P P'(-P)-(1+S)pg(-P)}, 

sin Its ^ 

-I < Rt Ci) < 0 


00 

{ (r^+d^-2drp)^'^^ 


jLi j^(.p) , 

sin us 


-i<r fhc5) < 1 

- 1 < p < '1 • 

The relations given below have been used to 
derive the above results. 


- ^-1 = ’ 

P P^{P) - 53.1^5 = S Pg(p) . 

(l-P®) = (s-t-l) p Eg(p) - 

(s+1) Vl^P) = + ® ’ 


p 


-m 

a 


(p) = (-1)“ 


fCs-m+l) Cp). 
T(s+m+l) 



(-D^d-p")"/" 


J- 

a(p)“ 


\ >* ... 

B(s, 1-s) *= SXIltlS * 

wliere B denotes the Beta fPnotion, m is a 
0 may he oofflplex and P = o°®9 ^7®'^ 9 < i> • 


positive 


integer, 



OMPTER 


vni 


miMERICJili EViiLUATIOlT MD DISCUSSION 


Ihe convergence of the integrals given in (119), 

will now te discussed. The asjrmptotic expansion of 

P ^ (cose) and P’ ^ (cose) may te derived from 
- 2 — 2 
the following results ((19)); 



P (s+m+l) 

T(s+ I ) 


£ 

{=0 


(-1)^ |-m) ^ 

\i (Sslne)^ (s+ !• 


. sin i (s+-{ + |)e + (| + |)» + |-t»] ( 162 ) 



1^9 


T (b+di) 

r(s+n) 


.m-n 


(m~n)(m+iv.l) p 

1 + + 0(s"'^) 

2s 


(153) 


as jsj -* oo, m, n real, -it < arg (s) < « . 

Using (152) and (153) ^following results may "be 
obtained . 


P . (p) 

~ ^ +ia 


ss 


e2cp(a.e) 

Y2ia~a 


{ Sig, + ] 


(154) 


Q S e S Ti - e f e>o, a** 


00 


and 

pi 


1 ' (p) » - 

- ^ +ia 


a 


exp(a.0) . 3p .2 


2% S7'^'' 


f 1 „ + 0(rt )}, (155) 


8ocq 


e4^0^Ti-C» G>0, a"* 


On tlie axis of tlie cone, i.e. for 0 *= 0 


P 1 (1) - 1 

- 5 '►ia 


and P’ (1) = -* I ^ i **■ 
- * +ia 


If each of the integrands in (119) is denoted by 
f(a, P» P )» then malcing use of the results (l54) and 
(166), it may be seen that 

f(a, P. Po) “ 0( ) 


as a 


( 166 ) 
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This coiapletos "biie proof for convergence* 

It seems necessary to find the stress distribution, 
for Some particular cases* Shis has been done by 
evaluating the integrals in ( 119 ) numerically for some 
cases. The method of numerical solution of these 
integrals may be divided into two parts. The first is 
the calculation of the Legendre fmictions of the type 

P . (p) and their derivatives. The series expansion 

- ^ + ia 

I* 1 (p) given in ((20)) and in (l57) below has been 

- §• + ia 

employed for calculation. 

^ 1 (p) = !’( ^ - ia;» ^ + ia ; li ), (157) 

~ + ia 

~1 < p < 1 . 

P* . (p) may be obtained directly. The series on 

- f + ia 

the right hand side of (l57) can be calculated upto any 
degree of accuracy. 

The second part is the evaluation of the infinite 
integrals. The presence of oscillatory functions sin ma 
and cos ma suggests the use of Pilon’ s method ((2l)), 
However, if the increment is small, one may employ even 
the well known Simpson’s rule. Pilon’ s formula is given 
below i 




b 

Let ss / f(a) sin ma da . 

^ a 

If tile interval (a, b) is divided into 2n intervals eacb. 
of length ■iy then, approximately, 


■i 


a^( f(a) cos(.ma)- f(b) cos (mb)) + 


where 

S„ 

2r 


03 

11 ]^f (a+2r'fv.) sin ^ m'(arf2r fv) j 


r=0 


- i |f (a) sin ma + f (b) sin 


S 


2r-r 


X” f I a+ (2r-l)i] sin m'J^a + (2r-l)'t^] , 


1 Sin 

a. =s — + 5 “• 

^ Px ^Pl 


2 


. 1 + COS 

®'53 2 

^ ^ pf 


sin 2p^ 


} 


( 158 ) 


^3 


4 ( 


sin 

_ 

^1 


cos 


) , 


where 


Tn' 


Simpson’s rule is rather Well known and is not 


given here. 



Ihe upper limit of tlie infinite integral may 
be taken as some finite value of a which is not 
arbitrary and is to be chosen suitably beoriug in mind 
the result in (156) * It may be remarked that value, 
of a may be suitably adjusted depending upon the 
semi-vertex angle of the cone, 

Ihe interval len^h is taken to be 0.1. 
lumericai work has been done for two values of the semi- 
vertex angles 0^ = 30^ and 45°. The values assigned to 

0 are 0° and q . The value of v) is taken as 0.25. 

0 

I'or each of the above cases the stresses , P^^ , 

Pcf<f and have been calculated. It may be noted 

that the resultant stresses are the sum of two systems 
of stresses ^^rr^2 first 

system is not complicated and may be calculated to a 
fair degree of accuracy. The second system of stresses 
involve infinite integrals and complicated integrands. 
Both Pilon's method and Simpson's rule, were employed 
and the results obtained by both the methods agree 
closely. It may be remarked here, that the boundary 
conditions are satisfied. 

The results of our problem for d =0.01 were 
calculated and compared with the well known results 
for d ss 0 given in ((16)) » It was felt that the two 
results cannot differ significantly. The values assigned 
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to r for d = 0.01 are 0.1, 1.0, 2.0, 3.0, and 4.0. 

It is interesting to note from latlo I, p.12>4 tliat tlie 
results agree closely witli maximum error of 1 percent. 

IMb provides a further check on the numerical solution 
of the problem. 

Graphs of , Pqq, surface 

as well as on the axis of the cones for = 30 and 45 
have been drawn. The values given to d are 1.0 and 8.0 and 
for each value of d, r takes the values 0.5, 1.5, 2.5, 

3,5, 4.5, 5.5, 6.5 and 7.5 . Pour graphs are given onp.l35y 
U'? nfl • are self explanatory. 
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OHAPSSR IX 

EI/ASEIO OYLiroER TJIDER M AXIiL POROE 

Tile problem of determining stresses due to a 
point force along the axis in an infinite elastic 
cylinder is considered in tMs chapter. The problem 
may be stated as follows : 

Let a point force (O, 0, 2^) act at the origin 
of coordinates in the direction of the axis of an 
infinite cylinder. The radius of the cylinder is a 
and its axis is taJcen along the z axis (-w < a < «>) 

(Pigure 1 ,p.155 )• The solution to the above problem 
has been achieved by the same process which has been 
described for the cone in chapter 6, p.loa . The first 
system of stresses are obtained by applying a point force 
in an infinite elastic medium. The solution of this is 
known from ((l6)). Note that in this case the point force 
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acts at (0,0,0). Ihe displacement components 
in Cartesian coordinates wlien a point force 2^ acts 
at (0,0,0) may le derived from (113) ly putting 
to be aero and (X^, Y^, 2^) = (0,0, 2^). Ibey are 


(A-f'/d) 2o 

1 , = ^ 

STCjJl- C A+S.}Jl) 

xz 

' ’ 

^3 

(A •+ 2o 

yz 

^ ” 8ii|i(A+a/x) 

3 ’ 

^3 


(159) 


u 


z 


CA+/^) 2^ z^ CA4B./x12^ 

d%ji( A + i yii) 


1 

^3 


, 2 „2 2 . 2 
where r^ - x + y + z . 

The displacement components in cylindrical polar 
coordinates (E, 9, 2) (x =* E cos0, y » E sin0, z = 2, 

0$ e < 2ii, -.« < 2 < “ ) are 

(M/L) 2^ E2 

"8 ’ ’ 

~ 8*/xa-*i/<-) (E^+Z®)^/^ 

Uq ss 0 . 


(E^+2^)®/^* 

(160) 
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Stresses corresponding to (l6o) in cylindrical 
polar coordinates are 

yu. Zq z nh 

(E^+2^)^/^ 4tiCA+2^^ 


(P 


/J. Z 


00 


” 47i(A+2/jl') (r2+Z^)V2 * 


R^Z 




5 Z^ 

4nCA4 3yL.) (R^+Z^)^/^ 


Z 


R" 


RZ 'l “ 471CA+3/1-) 


(3 A +4|i)Z^ 

4ti (A + a)^ (3+i^/2 ' 

(5A +4^) z, R 

4ii CA+iyuL) (r2+Z^)^/2 


CPj^)l « 0 and (P 02 ^i= 0 . 


(161) 


lElie autscript one is used to the above stresses 
in (l6l) to indicate that they are the first system of' 
stresses. 

Row if the point force Z^ acts at (0,0, O) in the 
cylinder, then on the surface of the cylinder, the normal 
and sheai^ing stresses are zero i.e. P^jj^ = 0, Pj^g" 

P =sOatR=a. So a second system of stresses 

(^br)^. C^EzJg- (^ee> 2 > <^22)3- <^ez >2 



superposed on the first stress system such that on 
the surface of the cylinder 


Putting R « a in (P^r)^, ^^RZ^l Cl6l), we get 



Che problem now is to find the stresses within the 
cylinder with the two boimdary conditions as given in 
(162} and (l63). To solve this, we start with the 
eq,uations of equilibrium in cylindrical coordinates 
(R,e, Z), Por axi symmetric case (variation with respect 
to e is zero) in the absence of body forces; they are 
given in (9l) und (92). The solution of these equilibrium 
equations in terms of two functions ^(R,Z) and -A (E,Z) 
is given in chapter V . The stresses in terms of (j) andJl 
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are given by (108) . ([) and Xi satisfy (112). 

Ilie exponential Fourier transform of a function 
f(R,Z), denoted by f CR,a) is defined as ((17)) 

f(R,a) = /” f(R,2) d2 (l64) 

—CO 

and tbe inverse Fourier transform is given by 


f(S,Z) = -i /” 7 (S.a) Ic . (165) 

2 % 

Assuming that 

gictZ _8_ and 8^“^ ((j), XL) vanish as a - ± “ 

az 

and tehlng the Eourler transform as defined in (164) of the 
sanations in (llZ), it may be seen that the solutions of 
the transformed eauations which are finite at E = 0, are 


-a a A(a) 


( 166 ) 


^ ^ { 2 B(a) I„(Ea)- A(a) H I- (E«)] , (16?) 

bar denotes the exponential Courier transform and dash 
stands for the differentiation with respect to E and. not 
with respect to the argument i A(a) and B(a) are the 
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coaistants to loe determined from tiie boimdary conditions, 
is the modified Bessel function of the first hind and 
order zero. 

Taking the Fourier transform of (^10^)2 
(108) and substituting the values of (j) and from 

(I66) and (167), it may be seen that 


2R(1-^J = A(a)[(2(l-0)+RV}i;(Ra)- 

+ 2B(a) [l'(Ea) - R I^CRa)] , 
2R(l--0) (^)2 * A(a) {(l-S'O )R a%(Ra) )Io(Ba)} 

- 2 B(a) I^CRa) » 

2R(l-^) )i'^zzh “ (2Ra^ a) + R^ 


+ 2 E E(a) a) » 

2E(l.^ )(^)g = «®l(a) I^CH oi) +2iHa B(a)Ii(Ea) , 

( 168 ) 

(i^)2 « 0 and ° * 

Taking the Fourier transfom of (l62) and (l6S), 


we get 





^ ^0 jua E^Calal) 

2% (A+5^) ^ , 

+ aoc KqC®'1cx1)|> 

(169) 

Jo 1 (A+57 ^)k’ aUl 

2% (A+4/^) L 

+ ( A+f)ii a^K^(alal)] ,(l70) 


K is the modified Bessel functiori of the second 
0 

hind and order zero . 

The following result has been used in (169) and 


(170) 



T(0 (2a)^“^/^ 


(171) 


given in ((22)) and recurrence relations in ((23)), ■ 

From first and fourth equations of Cl68) and 
the boundary conditions (l69) and (l70), the values of 
the constants A and B may be determined. They are 


A(a) 


2i 2, 


(l-^) )a r^(A-»-a.yd.) + a^a^CA+/^^} ^^^(alai )lQ(aa) 


- a^a^C A4-/A. XK^Calat) ^^(aa) + ( A+^/l ) 


(172) 



14 ? 


aiid 


B(a) 


i21o(l-.iJ) a 
7i;(A+a^)A (a) 


C (2(l-v^ )+a^a^j E^Cajai )Io(aa) 


-(A+f) aV lo(aa) Ko(alal) + ct\( A+ )+aV ( )! 


where A (a) = - 2aa 


Cl73) 

(2(l--v')+ a^a^} (l^Caa))^- a^a'^(lo(aa)fJ 

(m) 

Putting the values of A and B in (l68) and making 
use of the inversion formula (165), the second system of 
stresses come out to he 


(P™) = - -SV T /“K(«n^2-EV+IlV)Ii(Ea) 

BR'g 27C%1( A-+a/i-) o L 


. E„ Io(B«)j- B,(a) (l^CEa) - X„(Ea)l] ia. 


Zr 


CO 




- / 


e 9'2 2 *®E(A* 3 /^)° 




+ I- (Ba) B^(a) 


3±n a Z aa , 

' Ai (aT 




si 


^ s/C A*s/<- ) 


OQ 

/ j^A^Ca) (BIoCHa) + Ea Ii(Ba)l 


+ I, 


a da , 


,(Ea) Bj^^(a)] 





2 % 




/[e a Io(Ea)Aj_(cc) +Ii(Ea)B^(a)] 




wRere 


0 Z '2 


( 175 ) 


iLj_(oc) = - (( A4;^/^) + aVCA^yM^ )j a%(^a)Ii(aa) 

^ aV(A.4/l ) Ko(aa) I^Caa) + C A-+^^)a^, (176) 

B^(a) » - ( ) ( 2 (l-V) + a^a^i a%^Ca<x)l 2 (aa) 

-(X45|x. )aV loCaa) Ko(aa) + a^[( C A+/^ )J, 

(177) 

Aj(a) = .[2(1-'^)+ a^a^} a^a^l^oC^)} 

It may be noted that tRe stresses in tRe cylinder 


are 


■00 


zz 


^ER == ^^RR^l ■*■ ^^RR ^2 ’ 
^^00^1 ■*■ ^^ee^2 ’ 


^RZ 




■R 0 


0 and Pq2 = ^ • 


( 179 ) 
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On putting R = a in (l79), ^RE and come 
out to be ziero, which should be the case. Also as the 
radius of the cylinder tends to infinity i.e. when a 


each of the stresses ^^08^2* ^^23^2* ^^R2^2 

aero and the results of axial force in an infinite medium 

are obtained. 

It will now be shown that the integrands in (lV5), 
regarded as functions of a are free from singularities . 

Ihe singularities of the integrands may arise from the zeros 
of A I (a) and the poles of the numerators. It may be shown 
that the zeros of A| (a) do not exist for the range, of a 
from 0 to « . The argument is as follows : 


A, (a) =* 2(l--d) |I^(aa)j^+ (I^(aa)j^ - 




Using the relation 


> m+n °° j m+n+2i. 

I„(a) I„(a) = ( I ) I 


|i- |m4 -:-^ | n^'^ istS±^ 


2l 

( I ) ,(180) 


it may be seen that 

» ■ 

2 r-r 




(l^Cacc)} - (ID^li+l 


act 2^ 

■( -T ) 


finally, the e^jpression for (a) written 


as 
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A|(a) = 


2 4 

a a " 

<=o 


(l-<-20(l-~.) )-| ll.^ 

(l+^)(2+g )'*(ir) 




aa . 2^ 


( 181 ) 

The expression in the square brackets of (l8l) is 

positive for all integral values of so long as 0<'^< g . 

Therefore there does not. exist any real zero of the 

denominator except at a = 0. However, it may he noted 

that in the limit as a - 0; the integrands in (175) 

either tend to zero or to some finite limit but do not 
become infinite. 

The integrands in (l75) contain the functions 
K^Caa) and K^(aa) which are singular at a = 0. But these 
are multiplied by such functions that when a - 0 and the 
limit of the product is taken, it tends to zero or some 
finite limit. 


Making use of the following asymptotic expansions 


{(23)) 





' 1 

W ~ 1 

r 

-lb 

(4n^-l^)(4n^- 3*) 

K„(a) = 

-- e-“ 

S a ' 

1 + , 
li 

(8a) 

i2 (8a)^ 


( 182 ) 


and 



• » » 


oo 





4n®-l® (4n®-l®)(4n®-3®3 

ll (8a3 * li (8a)'2 



a , (183) 

the order of convergence of each of the integrands in 
(l75), may be given as 

0 ( — -)a ^ g^g Qj -* oo . (184) 


The integrals in (l76) have been evaluated nuraorically 
for some cases. The method of the numerical solution of 
the integrals given in (l75) may be devided into two 
parts. The first one is that of calculating the, modified 
Bessel functions lQ(a)» Ij^(a), K^(cx) and K^(a) . For 
calculating these functions, the series e^jpansions given 
below have been employed ((23)) . 


1 a n+2r 

I (a) a S” — — — ( - ) , (186) 

r(r+l) r(n+r+l) 


• ' f] . n-1 2m-n 

K^Ca) » 2 2 ^ 


OO 


msso Tf 

+ i(-l)^ S — — — — j ^Cn+m+l)+ I • 

(186) 


jm Ijm. 



For large values of the argument, asymptotic 
expansions for K^Ca) aud K^(a) given in (183) add (l84) 
have been used. She values of I^(a), I^(a), K^(cx) and 

calculated for different values of the argument were 
found to he the same as given in ((24)). 

She second part is the evaluation of the infinite 
integrals* She presence of oscillatory functions sin aZ 
and cos aZ , suggests the use of L.N.G. Filon’s method . 
However as stated in chapter Till, if the interval length 
is small, one can employ even Simpson’s, rule . 

She upper limit of the infinite integral may he taken 
as some finite value of a which is to he chosen suitably 
keeping in mind the result given in (l84). She infinite 
integrals were truncated at a = 20 and « = 25 and the 
difference in the values of the integrals .for these two 
truncation limits, is insignificant. She interval length 
was taken to he 0.1 . She ratio z/a takes values 0*0, 
l.Q, 2.0, 3.0, 4.0 and 5.0 . She value of Poisson ratio 

is taken as 0.25 . 

It may he noted that the stresses in the cylinder 
are the sum of two systems of stresses “, the first and 
the second stress systems as given in (l79). She first 
system of stresses are not complicated and can he 
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calculated witli a fair degree of accuracy. The second 

system of stresses involve infinite integrals. Filon* s 

method and Simpson's rule "both were employed for this 

finite interval and the results obtained by both of these 

methods agree closely. 

It is interesting to note the behaviour of 

^ee ^ZZ z . it may be seen that 

^RZ ^00 Z -* <» . But P 22 ^ 

constant value as Z -♦ <» . The integral for P^g 

form . ’ 

<90 Sin aZ ■ 

/ f(a) — da 9 
0 a 

which has the value | f(0) as Z - « . The analytical 
value of P 22 is -Z^ /2» a® as Z - ■» . Ihe nuaerioal 
value of P^z is -o■ 4 ')'^ 8 1 » ZoAa^ whioh Is attained 
approximately for all Z/a > 4 . 

The graphs of stresses Pj^, P|^2 » ^00 ^ZZ 
cylinder have been drawn for all t:^ cases mentioned 
above and are given in figures Z, '4,^ and S in the 
Appendix following this chapter. A comparison of these 
stresses has been done with the stresses due to a point 
force along the Z axis at the origin in the infinite 

medium in figures 3, 5» 1 ^ • 



iPPEiiDIX TO CHAPTER 


IX 




156 



157 


-I 









158 
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160 


CEAPTEE X 

POINT FORGE II m IIPIIITE ELASTIC WEDGE 

In tills chapter the following plane strain 
problem is considered. 

Suppose there is an infinite elastic solid, 
which is bounded by the two planes 0 = + 9q (R,e,2 
being the cylindrical polar coordinates) . Thus the 
edge of the wedge is along Z axis. In Cartesian 
ooordimtes, therefore, the x, axis is the axis of 
symmetry and the ^ axis is perpendicular to Z and x 
axes. Iota that i axis in Cartesian coordinates is 
coincident with Z axis in cylindrical polar coordinates. 
The system is supposed to be in a state of plain strain. 
Hence we might consider only a plane cross-section 
perpendicular to the Z axis. We shall henceforth work 
in Cx,y) pl««ie only and consider the stresses ^^x’^xy 



mid or tlieir equivalent in polar coordinates. 

It may te remarked at this stage, that the 
probloms worked in the previous chapters wore axi- 
aymnietric , But in the present problem the point force 
is not necessarily acting along the axis of symmetry . 

Let a point force (X^, Y^) act in the wedge at a point 
whose Cartesian coordinates are (h,k) (Figure 2 p.l39). 

Iho solution to this problem has been acliieved by the 
superposition of two stress systems. The first system 
of etresseo is that due to a point force in an infinite 
two-dimensional elastic plane. This we have named as the 
first stress system. On the faces of the wedge 0 — + 0 q> 
normal and shearing stresses may be calculated from the 

first stress system. 

me second system of stresses are oMainea ty 
BOlvlnS the followins proUom : consider a wedge whose 
faces are hounded hy the lines 6 = i (note that we are 
worhinS a two-dimensional problem) . -ipply '*'he 
of the wedge end which are equal in 

magnitude hut opposite in sisr to those obtained from 

first stress system at 9 = + Sp • 

to the wedge. These stresses are the second system of 
strousos. The superposition of these two stress systems 
would give the stresses due to a point force rnsxde the 



wedge, v/ith zero tractions on the boundary. 

Ihe first system of stresses i.e, when the point 

force (X^, Y^) acts at a point (h,h:) are given by (llA) . 

aihese stresses in polar coordinates (r,e) (x = r cos e, 

y a r sin e) may be easily derived from (114)- by 

transforming the stress system in polar coordinates. 

They are 





4ic(l-.v) 


(1-2^) 


£r cose - h cos28 -ksin2ej 



-»■ 2(r cose-h) -(.r-hcose - k sine] 





j 

Y 

0 

r fi p^) ^ ^ cos 2 ej 

1 

I 


2 

^ 2(r sine ~k) {.r- h cose - k sine} 



X r ' 

0 1 

(l-2v)) {,-r COS0 + h cos2e + k sin20] 

^ iL ^ — - * 

4Tt;(l--vi ) ^ 


2 

2(r cose- h) (k cose - h sine) -j 


t>4 -* 

Y 

0 

■ (i-2^) ^-r sine + bi Bin26 - k‘ cos 20} 
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2(r sin0 - k) (k cos0 - h sin 0 ) 




] 




X r (l-S'^)) ir- ^ sin0 + k sin26 - k cob20^ 

0 1 g- 




r: 


2 r'^coso Ck siiie - k cosej^ 


r: 


2r (2hk cos^e + sineCk^cos^0 - k^- k^cos 0 )j. 


R. 


^kk cos’20 + sin0 cose(k^~k^)} 


2k 


R. 


Y f (l— 2'vi ) COS0 — k cos20-ksiii20j 


4ti:U“V ) 


irrr [ 


Et 


2 r^sine (k sin 0 - k cos 0 ) 


R. 


2r I cos0(k^+k^sin^0 - k^sin 0 ) -» 23il£ sin© } 


R. 


2 2 ^ 

(kk cos20 + sin© cos0(k -k )j 


2k 


r: 


(187) 



where + h^ + - 2r(h cosQ + k siro) , 

The subscript 1 refers to the first system of 
stresses. Prom the above expressions of .stresses given 
in (187) , and at e = + may be easily 

calculated^ Rovir apply on the faces of the wedge normal 
and shearing stresses and P , equal in magnitude 

and opposite in sign to those obtained from (187) by 
putting 0 K + 00 . 

Thus, v/e apply (Pqq)^ and (Pj^q)^ on the faces of 
the wedge such that 


{^^60^2 ' 

J 0 =: + Bq “ " -[.^^06^1^0= + 0Q 

(188) 


: e = + 00 " “ e = ± ©o • 

(189) 


The solution of this auxiliary problem will now 
be described. 

The equations of equilibrium in polar coordinates 
in the absence of body forces and their solution in terms 
of a function ^(r,e) which satisfies the biharmonic 
equation 



\ — ^ 

ar r 



>2 

.2 


a~ a'^t) 

)( 


r^ ae^ '' ar^ 


1 

r 


d([) 1 8^(j) 

+ ~~o '* 9 ' ) = 0 

ar r ae'^ 


( 190 ) 



is given in details in ((17)). It is sufficient to 
mention here tliat the stresses may he expressed 
directly in terms of the function ({) such that 


1 ii ^ 1 ^6 

■" + p "o' » 

r 0r T 00^ 








(191) 


Now if ^ is defined as 

— OO 

^) = / ij) dr , (192) 

0 

then it may he seen that the solution of (190) is 

(|) *= A sin(s-2)0 + B cos(s--2)e + 0 sin s0 + B cos s 0. 

(193) 

This solntion in (l93) may also he obtained from 
the solution given in ((17)) hy putting s~2 for s (s 
has been replaced by s-2 for convenience). 

Multiplying each of the e(juatioxis in (l9l) by 
and integrating with respect to r from 0 to and 
using (193), we get 



1 (i 8 


(P = - A(s-2)(s-l)sin(s-2)e - B(s~ 2 ) (s-l)cos(s-2)0 
- C(s+2)(s-l) sin B 0 - D(s+ 2 )(s-l) cos s 0 , 

(Pqq)^ = (s-.2)(s-l)( A Sin(s-2)e + B cos(s-2)e + 0 sin s0 

+ I) cos s0] , 

(P^q)^ = (s-2}(s~l)^ A cos(s-2)0 - B sin(B-2)0j 

+ s (s-l) {c cos s 0 - B sin sej , ( 194 .) 

wliere bar denotes the Mellin transform of a function as 
defined in (l26) and the subscript two has been used 
for the stresses in (194) as this will be the second 
stress system. 

Multiplying the boundary conditions in (188) and 
(l89) by r®*"^ and integrating with respect to r from 0 
to , we get 



0= + 00 



dr 

Q* t ©0 



( 195 ) 


and 
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0 =+ 0 , 


- / r 
o 




0= + e 


dr 


Or 


(.^^xehl 0^ + 0 


‘0= + 0 


Tiie infinite integrals occurring on the ri^t 
hand sides of (l95) and (196) may he evaluated using the 
result 


f -u.g-.£ — ™ — dr — __ = - % a®“^cosec 0 cosec ( tis)* 

0 (r +2ar cos0 + a) 

. sin (s-.l)e (197) 

0< Re(s) < 2 , a > 0, 

-u < 0 < It , 


given in ((l8)), after suitable modifications. 

bet the right hand sides of (195) and (196) he 
denoted hy respectively, so that 




0 = 0 , 


fg2(s,eo) > 




0 = -0, 


f22(s> -eo), 




0 =i 0, 
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^ “00^ • (198) 

0 = -Gq 

Tile ei^pressionisof f^^(s, e^) , f^^(s,-0Q), ^ggCsje^) and 
fgg^®’ “Qq^ real and imaginary parts are given in the 
Appendix II following this chapter . 

It may be seen that all the integrals on the 
right hand sides in (l95) and (196) are valid for 0< Re(s)< 1. 
With the help of (l94) and (l98) , we can determine the 
four constants A, B, 0 and B. They are 

|f2_^(s,eQ)+ f^^(s,-.9Q)| (s-.2)sinseo - (fggCsjeQ) 

^ __ -fggCs ^-epDjs cos S9 p ^ 

S(s-2}(s~l) {(b-I) sin 20^- sin 2(s-l)0pJ 

(199) 

|f^^(s,-ep)- f^^(s,9^)j (s-2) cos SG - 

+ ^ggCsj-Gp) }s sin sGq 

B S3 - "■ ' ■ ■ ■ '■■■ ■ ' ’ 

2 (b- 2 )(s- 1}{ (l-s)sin 20^- sin 2(s-l)e^J 

(200) 

jf^^(s,9o)+ f^i(s,-0g)}sin(s-.2)0^- 

C =s - ^ ^ ^ ’ 

2(s-l} {^(s-l),sin 20^ - sin 2(s-l) 63] 


(201) 



D - 


_|f j-Gq) — f s )0 q)| COS(y~2)0Q~ ,0 q) 

+fpp(s,-.eo)jBin(B-2)eQ 

. i .(202) 

2(s-l) ^ (l-s) sin 20^ - sin 2(s-.l) 0^^ 

Substituting tliese values of the constants 
A,B,0 and I) in (194) and inverting (i94) vdth the help of 

given in (l27), we get 

r ^(s-2) sin seQsin(s-2)e 

L A,(s,0o) 

^2 f 

S0 sin(s-.2)6_\ - is cos sO^ 

4i(s,eo) 

» 

0. 

.sin(s-2)0 -(s+2) sin s0 cos(s-2)g 1+ 

A2.(s,.0o) 

. ((s~2) cos S0Q cos( 3~2)0 » (s+2)cos S0 cosCs-2)0Qj 

I 

^2 ( s sin s0Q cosCs-2)0 - (s+2)co3 00 

j0q} 

sin(s-2) ^olj * 

(V ) « J- / ^^(s-2)(sin se sin( 0-2)0 - 

cnee’s = 

- sih(s-2)9Q sin s0] - 


tlie inversion formula 




rr''2 


y+i" 

i / 

27ii ■v'-i'=^ 


-(s+2) sin 





wli^re' 

(s,e 


Og s COB s0Qsin(s-S)0 - (s-2)cos(s-2)0oSi3is0] 

^i(b,0q) 

(s-2)^cos s0qCos(s~ 2)0 - cos s0 cos(b-2}0q^ 


4j^(b,0o) 


Gg _^s sin s0qCos(s-2)0 -.(s-2)cos s0sin(s~2)eQ| 


ot£. 


A 2 . (s»0o) 


1 C 


(203) 


StcT 


f 

I 


1 I ''o 


2)sins e cos(3-2)0 


V-i°° L A , (sjBo) 


0. 


, •, '2 s(cos s0^cos(s-2)0 

- s cos 80 sin(s-2)e j - ;; . G 

A I (s,0o) 


* 

, , ;-Cs-2)cos s0„sin.(s-2)0 

- cos s0 cos(s-2)0 J + — i 

A S r 0 o ) 


+ s sin s0 008(8-2)0^1^+ . 

vS J 00^ 


1 

2 s|sin s0QSiii(s-2)0 


sin 80 sin(s~2)0oj 


ds , 


s 2 |s sin 20q - 2 sin 80^ cos(s-2)0q|, 


(204r) 




Ai(£j,0o) = 2 s sin 20^ - 2 cos sOq sin(s~2)eQ|, (205) 


" f^^(s,-0Q)- f^^(s,0Q) , 

( 206 ) 

Cg — fgg(sj0Q) — ^^^(sj—Bq) ; 

*^2 ^ 22 ^®’® o ^ ■*■ • 

Th.e integrals in (203) are not affected by the 
particular choice of V, as long as Re(s) (= 7 ) remains 
within the limits 0 < Re(s) < 1 . The choice of V = l/2, 
which in turn iic^lies 

sr= i+ia, , (207) 

is suggested Toy the fact that the functions occurring in 
the integrands of (203) are simplified to a great extent. 

In order to establish thfef -validity of the choice 
of V , we have yet to show that the integrands in (203), 
regarded as functions of s, are free from singularities, 
if (207) holds. The singularities of the integrands may 
8 , rise from the zeros of sin sn which occur in the 

denominators of f^^Cs,©^), f^^(s,- 0 Q), f 22 ^®»®o^ ' 

fgg(s,-. 0 ^) and the zeros of Aj(s, 0 q) and , It 



way be noted that sin sit does not vanisli for the values 
of 0 ^,lven in (207). A, (0,60) A:^(G,eQ) may be 

written as 

( 0 , 0 ^) =: - si3:j2eQ+ 2sin8QCOsh 2aeor-2i(coseQSinh 2a0Q 

- a sin 200 ) , 

(s,0o) = - (sin 200+ ain0ocosh 2 q: 0 o)+ 2i(a sin20o+cos0o, 

sinh 2a Bq) . 

It may be seen from the above eixpressions that 
Ai(b, 0 o) and Aj^(s, 0 q) are not zero for positive real 
values of a and 0o > 0 . 

We nov/ put s = i. + ia in (203) and separate the 

2 

integrals in real and imaginary parts. The real and 
imaginary parts are found to he even and odd functions of a, 
respectively. Hence, the contribution of imaginary part is 
zero and we obtain the complementary stresses as given below 




JJ3 


I 7 § f/W~ ^ “^1“ “ ^22®^^ ^^a)d(x , 


%T ' r 




00^2 - 'J'A 


00 

j (fggcos m^a ~ ^44sin m^a) da » 


%x'"' "Tq ' ^ 0 




ggcos m^a - I'ggSin m^a) da , 


( 208 ) 



where 


= log (ro/r ) , = h® + ir®. i 

the (luantities S'gg, and are giyen 

in the Appendix II to this chapter. 

Pinally, the stresses in the wedge due to the 
point force are given by 

^rr * ^^rr^8 ’ 

^00 " ^^00^1 ■*' ^^00^2 ’ (209) 

^r0 “ ^^r0^1 ^^re^2 * 

where etc. and iS etc. are given by (187) and 

(208) respectively. 

It may be verified that on the faces of the wedge 
0 a + 0Q , the normal stress P^^ and shearing stress P^^ 
vanish, as they should. On putting (209) , 

the results due to a point force in the half plane are 
obtained. It may be remarked tlmt while finding the 
second system of streesee In the wedge another problem of 
■that of preeorlhed tractions on the faces of the wedge is 

also solved. 

The IntagralB in (208) have been eTrel-uated 
numerically for some oases and the numerical evaluation 
and discussion is given in the Appendix I to this chapter 
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APPEEDIX I TO CHAPTER X 

The integrals in (208) are all convergent and if 
each of the integrands in (208) is denoted hy f(a, 0 , 0 Q), 
then 

f(ct, 9, e^) = 0 ( a (e“(e-9o-ei) + e“(6-8o-92^ )}, 

h cosGo+k sineQ k coseo-ksin 0 Q 

where cos 9- = - — ^ 5“177>‘'' 6 p = - p-— p* ^ 7 - 0 — . 

The integrals in (208) have "been evaluated 
numerically for some cases. The infinite integrals may ije 

truncated at some finite value of a, hearing in mind the 
order of convergence given above which depends upon a> the 
point of application of the force and the semiangle of the 
wedge. The semiangle of the wedge 9^ has been taken as it/S 
and 7 c /4 and the point of application of the force (Xq,Yq) 
is taken as (l,0). Two cases, when the point force is 
acting along the axis of symmetry and when it is acting 
perpendicular to the axis, have been considered. The upper 
limit of a is taken as 25. Poisson ratio is taken as 0.25. 
The methods employed for the numerical evaluation of the 
Integrals are the same as described in previous chapters. 

Boundary conditions were satisfied numerically 
also. It was felt that if the point of application of the 



force is very near the vertex of tlie wedge, the stresses 
should not differ substantially from the stresses when 
the point force acts at the vertex of the wedge. Ihe 
solution of the problem when a point force acts at the 
vertex of the wedge is given in Love's book ((16)). 

Stresses were calculated for our problem when a point force 
along the axis acts at a point 0.01 on the axis of the 
wedge of semiangle 7c/4. These stresses were compared with 
the stresses due to an axial point force acting at the 
vertex of the wedge of semiangle x/l which can be calculated 
from the results given in ((l6)). This comparison of the 

Orf. 

two systems of stresses so obtained for some valuesir and e 
is shovm in Table 1 p. ni . As expected, the stresses do 
not differ significantly and the maximum difference is 
This gives a further check on the calculations, numerical 
values of stresses Pqq and when a point force 

acts at (l,0) perpendicular to the x axis in a half plane 
are given in Table 2 , p.lva > r varies from 0.5 to .7.6 
with an increment of 1.0 and the values given to 0 are 0 q 
and 0 . 

Variation of stresses P^^, P^q and P^ is shown 
in the form of graphs from 1 through 6 which are given in 
this Appendix from p . 119 to 184 . The semiangle of the 

wedge is 7i/4 and the point of application of the force is 
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I'ilBLB 1 

POUT POROE = (X^, O) . 

POIilT OP iPPLIOATIOlI OP TBE PORCE = (O.Ol, O) . 

= u/4, ; 'J = 0.25 . 

P „ aiid P „ are zero each. 

00 re 


0 

r 

1 

" IrT pr^erit 
problem. 

I* /K 

r^r, 

Ill solution 

0 

0 .1 
0.2 
0.3 
0.4 

.-7.7875 
-3.8893 
-2 .5930 
-1.9448 

-7.7797 

-3.8898 

-2.5932 

-1.9449 

%/Q 

0.1 

0.2 

0.3 

0.4 

-7.1868 

-3.5936 

-2.3958 

-1.7969 

-7.1875 
-3 . 5937 
-2.3958 
-1.7969 

%/4; 

0 .1 
0.2 
0.3 

0 ,4 

-5 .5040 
-2.7517 
-1.8340 
-1.3754 

-5.5011 
-2.7505 
-1 .8337 
-1.3753 




TABLE 2 

POUT PORGE =: (0, ) 

POINT OP APPLICATION OP THE PORCE = (l,0) 
0Q = tc/2 , V = 0.25 

P^„ is zero identically. 

00 


0 r 

^ An 

P /I 
r0^ 0 

0.5 

0 . 0000 

0.2963 

1 .5 

0.0000 

—1 .3440 

2.5 

0 .0000 

-0.5053 

= 0 3.5 

0.0000 

-0.2919 

4.5 

0.0000 

-0.1932 

5.5 

0 .0000 

-0.1379 

6.5 

0.0000 

-0.1036 

7.5 

0.0000 

-0.0807 

0.5 

.0.1600 

0.0000 

1.5 

-1.6331 

0.0000 

2.5 

-1 .6886 

0.0000 

3.5 

-1.4254 

0 .0000 

s: 71/2 4,5 

-1.1909 

0 .0000 

5.5 

-1.0109 

0 .0000 

6.5 

-0.8739 

0.0000 

7.5 

-0.7677 

0.0000 
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®11 = - I sinCe+Go) sin ( cosh aCe^^- O) 

0+0 

+ I sinC-^°) sin (0^-0) cosh a (0^+0) 

0+300 0 - 30 ^ 

-2 cos( — g — ) cosh aC0Q-0) + 2cos(--2 — )cosh a(0+% ) 

©o**®'* . / ^ 
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from the expression of f^^(a) 0 q) by replacing 
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i Y 
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0 
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the expression of 

,0i.by02 %2 
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fggCa. -60^ obtained from 

fggCa, Bo') replacing e„ by -6^ 
(a, Bo) *11’ *Z2’ 



CHAPTER XI 


CIRCULAR INCLUSION IN A VraUGB 

The method described in chapter X can also he used 
to solve soma other interesting prohlems , This has been 
done in this chapter by solving the problem of circular 
inclusion in a wedge. Ihe problem may be stated as follows: 

Consider the infinite elastic wedge described in 
chapter X. The system is in a state of plain strain. Thus 
we are considering only a plane cross-section perpendicular 
to the z axis of an infinite elastic solid, which is bounded 
by the two planes 9 = + z (=? Z)being the cylindrical 

polar coordinates) ; the edge of the v;edge is along z axis. 
The X axis is along the axis of symmetry and the y axis is 
perpendicular to x axis^ We shall henceforth work in (x,y) 
plane only.' The origin is taken at the vertex of the wedge 


and is named as 0 . 
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Let there be a circular hole of radius r’ in this 
infinite wedge (Figure 3, p.130 ). 3?he Cartesian components 
of the centre of this circular hole are (xq, . Suppose 
a circular disc of dimensions slightly bigger than those of 
the hole is embedded in this hole. I>ue to misfit, stresses 
would develop everywhere in the wedge. The problem considered 
is that of determining this elastic field so developed. 

The solution has been achieved in the following 

manner . 

Consider a circular inclusion of radius r’ in a two- 
dimensional infinite elastic solid. The centre of circular 
inclusion is at a point 2 ^ « + i y^ with respect to the 

origin at 0. The boundary of inclusion is denoted by L . 

The inclusion in the absence of the surrounding material 
undergoes a displacement whose Cartesian components are 
characterised by (G^x, G^y) with respect to the origin at 
Zq ; the centre of L . 

Let ( u"^ , v*^) be the displacement components in 
Cartesian coordinates of the inclusion and (u , v ^ be 
those of the matrix, where + and - superscripts have the 
same meaning as in chapter I. At the equilibrium 
interface if the displacement components for inclusion 
and matrix are denoted by , v^'*’) and (u^ , v^ ), 

respectively , then 



^02 


- "b~ = - . 

- V^' = - Bj^y ■ (210) 


Ihe sectionally xiolomorphic functions <j)^(z) and 4^^ (z) 

(z = X + iy = r-e,^®) having the line of discontinuity L 
may be found out as discussed in chapter I. It may be 
noted that the sectionally holomorphic functions 
and (z) are determined with respect to the origin at 
Zq . The sectionally holomoiphic functions ^(z) and 
'4' (z) with respect to the origin at 0 may be obtained 
from f)^(z) and (z) by the translation of origin ((-2)) 
and are given below. 


Mz) == (|)^(z-Zq) , 

H'(z) = S^,(z-Zo) - Zq‘ (|)^(z-Zo) . 


( 211 ) 


Stresses and P^^ in polar coordinates 

(r,0) can be easily calculated from (21l) using (lO). 
Let these, stressed be called the stresses due to the 
first stress system. They are denoted by (P ) , 
^^06^1 subscript 1 is added to the 

stresses to indicate that they are the first system of 
stresses, calculated at 

0 = + 00 • Now apply (Pqq)^ and (P^}^ eciual in 



magnitude and opposite in sign to (^09)^ 
tile faces 6 = + Gq of tiie wedge. Tims 



This will develop stresses everywhere in the wedge . 
These stresses are called the second system of stresses 
and we distinguish them hy adding subscript 2 to them. 
Superposition of these two systems of stresses ihside the 
wedge would give the solution to the problem when the 
circular inclusion which was hitherto present in the 
infinite plane, is situated in an infinite elastic wedge 


bounded by the lines 0 = ± Gq * 

We now deteimiine the first system of stresses. When 
a circular inclusion of radius r’ whose centre is at a 
point Zq Czo « Xq + iy^) is present in a two-dimensional 
infinite elastic plane, the sectionally holomorphic 
fimctions (j)(z) and ^ Cz) for inclusion and matrix are 

respectively 



2/lG^ 


(z- 




9 


(K+1) 



+ 


2 p. e ^ 
(K+l) 


(212) 




z 


0 


and 


ll)„(z) = 0 




4 fJ. e^r 
(K+l) 


.8 


(z-aj 


(SIS) 


where the su'bscript i and m refer to inclusion and matrix 

respectively. It may be noted that t)j_(z), 'V^Cz), 

and are determined with respect to the origin at 0. 

Ihe stresses in polar coordinates for matrj.x and 
incD-usion corresponding to ^(z) and 4* (213) and 

(21 2) are 




(K+l) 


4^e^r 

(K+l) 


f2 


/ 2 2 
|r +r 


11^°^ 
^r^+r*^ 


2(e“^^l)--2rriiC0 3 


11- Srrj^^oOBCe-p^j^)} 





4>J.e^r'® Bin 2(6-p :|^)-araj^iSi=i(9 -^ii)S 

(K+l) oos(e-Pj^]^)j® 


and 


( 214 ) 



aos 




4 ^ 
(K+l) 


Fee^iito (K+i) ’ ^ H 


"1 

where (^0/^0) 


(215) 


On the faces of the wedge 




4)i.G^r‘^(r^+r^^cos 2(0Q~^^^)-2rr^^cos(6o-£^j 


0=00 (K+1) 2rr^^cos(eo-^^l)j 


r 


i^^66^ale=_9, 


( 216 ) 

4 /JLe^r*^ Ir^+r^l cos 2(eo+^^j_}“2rr^^s(eo^^ 


(K+l) 


2rr^^cos(0o4p^^)j 


(217) 




0 = 0 , 


4fLe^r’^ 2(0o-p^- 2rr^j^sin (0o-. pj_^)j 

(K+l) {r^'+ - 2rr^^cos(0o-?il)j^ 


(218) 


•l^^r0^2\ 


0 =- 0 , 


4 fJLG^r»^[-r^^sin 2(0 ^+ p^ ^)+2rr^^8in(0o+^^ }j^ 

(K+l) ('r^+r|^ - 2rr^^ cos(0o+Pil)j^ 


( 219 ) 
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The solution of tliio au:d.lio,ry problem of a wedge 
with the boundary conditions given in C2l6)-(219) can be 
obtained with the help of the method given in chapter Z. 

As mehtioned there, we start with the equations of 
equilibrium in polar coordinates. Their solution in terms 
of a biharmonic function |)(r,0) is given in details in 
((17)). Stresses in terms of (}) are given by (l9l). If f 
is defined as in (192), then solution of (190) is given by 
(193). Also (P^)^, (P^q )2 and (P^q) 2 , the Mellin transforms 
of P , P^cv and P _ in (191) are given by (194). The 
constants A, B, G and D in (l94) are deterroined with the 
help of the boundary conditions (2l6) - (219). Multiplying 
the equations (216), (21?), (218) and (219) by r “ and 
integrating vdth respect to r from 0 to ■» , we get 

0=*y 0 

(220) 

-(.(WsL 

■^y= *““0 

f 22 (s»’-eo)» fiq(B,eo) and 


where 
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rofil and iraa^jinary parts are gi-ven in {22M) "belov/. 

It may be noted that all. tlie integrals on the 
right hand side of ( 220 ) are found to be valid for 
0 < He (s) < 1 . 

The four constants A, B, C and D may be determined 
from (194) and (220) . They are given by (199), (200), 

( 201 ) and ( 202 ) where fii(s,eo)» f^^(B,-eQ), ^gCsjeo) 
and f 22 ( 0 ,- 60 ) different for this problem. Substituting 
these values of the constants A, B, C and I) in (l94) and 
inverting (194) with the help of the inversion formula 
given in (127), we get (203). It may be noted that the 
values of 0 ^^, , Gg and Gg are different for this 

problem as ^o^’ f 2 ^^(s,- 9 o)» ^ 22 ^®’ 

are different for this problem. Other quantities remain 

the same . 

The choice of V — l/2 which in turn implies 

a = I + ia 
2 

is justified by the same arguments as given in chapter X, 
p. I,tl 4 . Putting s = I + ia in (203) and separating them 
in real and imaginary parts, we get the second system of 

stresses (^ 37 ^ 2 » ^^ 99^2 ^^r 9^2 





J (P^^cos mga 
o 


Pggsin mga) da , 



^^ 00^2 




(S'ggcos Biga - ^'440^^ raga)< 3 .a 


^^re^2 


1 




j>°“ C^55 “-2°^ “ ^66^^^ ni2a)<3.a » 

o 

(221) 


2 2 2 

wliere mg *= log(r4^/r) , , 

the quantities I'gg, S'44, ^'55. ^66 

as given in Appendix II to chapter X . 

!Pinally , the stresses in the wedge for piatrix and 

inclusion are given ty 


« 



^^rr ^2 ’ 




1 



( 222 ) 

tt 

•H 

U^rr^il^ '*' ^^rr ^2 » 


II 

•H 

CD 

CD 






( 223 ) 





The exproGGioiifj of f^ 3 _(s, 9^), 
f 2 g(c, 0 q) and ^^Cs-j-Gq) for tMs problem in real and 
Imaginary parts are given below. 




2 JULTI 


*2 


3 

(K+l)sin ©j^cosb. %a 


2 si^(9o-Pil^ { 


1 

Sin ^•^sin e^. cosh a ( 71 - 8 ^) 

0 - Q-i 

- a sin©^ cos ^^sinh aC7t-0^) - cos ^ cosh aCTc-e^)]- 

1 30^ 

+ sin 2(0 q-P^^)| ^ sin G^sin ^ cosh a(7C- 9^^) 

30J 9rf \ 

+ a sin e.cos sinh a(Ti:-e.)+ cos ^ cosh a(7C--0^)J 


where 


Zi^% r 2 sin(e^-^^^) ^ - sin 

y . \ . ^ ^ 1 . . . “ 


(K:+l)sin^0^cosh %a 

01 

sihh aC^t-Qj^) + <x sin ©i^sin ^ cosh aC 7 i- 0 i) 

0 

+ ^ sin 01 cos sihja ai%- Qi)]"*' 

0 S©! 

j- sin sinh aU-ei) - a sin ©isin cosh aCTc-Gi) 

+ I sin ©1 cos sinh aCTi-ei)^ J » 


cos ©1 = 


XqCosGo + YoSinGo 



(xq + yo)>- 
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-Oq) may be obtained from the expression of 
f^^(a, Gq) by replacing Gq by -Oq and by 6g where 


cos Q2 = 


Xq cosGq - To 

( E 2^l/2 

+ y,) ' 






,|2 


(K+l)sin'^0-| coshTta 


0 

!• sin ^7^cosh a ( -rt-G ^ ) sire - 


0 30 ^ 

+ a sine^cos g^sinh a(it-ej_) + cos ^ cosh a(7c--0jL) 

0 

+ 2 cos(0o-^^i)| I sin g-^sine^cosh 

0 0 -l 

- a sine^cos ^^sihh a(ii-B]_)- cos g cosh a(7t-e]_)j 

30 

+ co0S(eo-p3_l){ I sine^sin -g-^cosh a(ic-0^) 

30 ' ^ "1 * / 

+ a sin ej_oos sinn o(ii-e^)+ cob g- cosh aCn-e^)] 

2 ®G-i 

2iii|i,e^r’ ^ j^sin ‘ g^sihh a(Tt-e^)+ asine^' 

(K+l)3i]?0]_cosh Ucc 

Q 

aln ^ oosh - !■ slnDj^oos g-^slnli aCit-e^^) 

+ a 003(00-^ I - sin l^sinh +■ 
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0 - ^ ^ ^ 

+ a sin e^sin-^ cosli aCii-e^) + ^ sin Q^sos ^ 

/ ^ 

giiih, a(Ti- 0 ^) + oos - sin ^ sinli 

301 ^ 1 

- a sin ej_sin ^^cosh a( 7 t- 0 ^) + 2 61°°® 2 

siiih. aCn-OjL) (2*2.4) 

f . (a -0 ) insy olstained from tlie e^^ression of 

fgg(a, 9o) ty replacing e^, by -Bq an^ 6i ®2 

It may be verified from C 2 W) ai^ (22i) 

_ _ ^ a — + 0 , normal and sliearing 

the faces of the wedge 0 - i Wo> 

Streesee vanleb as they ahould. Ihe continuity of normal 
and shearing stresses at the inclusion boundary may also 
be verified. The results of circular inclusion in a half 
plane may be obtained by tahing the semi angle of the wedge 

as 1 !c/2. 


/ 
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